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In a variety of models the motion of the inflaton may trigger the production of some non-infiaton 
particles during inflation, for example via parametric resonance or a phase transition. Such models 
have attracted interest recently for a variety of reasons, including the possibility of slowing the motion 
of the inflaton on a steep potential. In this review we show that interactions between the produced 
particles and the inflaton condensate can lead to a qualitatively new mechanism for generating 
cosmological fluctuations from inflation. We illustrate this effect using a simple prototype model 
g 2 (4> — 4>o) 2 X 2 f° r the interaction between the inflaton, <j>, and iso-inflaton, x- Such interactions are 
quite natural in a variety of inflation models from supersymmetry and string theory. Using both 
lattice field theory simulations and analytical calculations, we study the quantum production of x 
particles and their subsequent rescatterings off the condensate <j>(t), which generates bremsstrahlung 
radiation of light inflaton fluctuations 8<f>. This mechanism leads to observable features in the 
primordial power spectrum. We derive observational constraints on such features and discuss their 
implications for popular models of inflation, including brane/axion monodromy. Inflationary particle 
production also leads to a very novel kind of nongaussian signature which may be observable in future 
missions. We argue that this mechanism provides a simple and well-motivated option to generate 
large nongaussianity, without fine-tuning the inflationary trajectory or appealing to re-summation 
of an infinite series of high dimension operators. 

This work is dedicated to the memory of Lev Kofman. 
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I. INTRODUCTION 



In recent years the inflationary paradigm has become a cornerstone of modern cosmology. In the simplest scenario 
the observed cosmological perturbations are seeded by the quantum vacuum fluctuations of the inflaton field [l[ . This 
mechanism predicts a nearly scale invariant spectrum of adiabatic primordial fluctuations, consistent with recent 
observational data [2j . In addition to this standard mechanism, there are also several alternatives for generating cos- 
mological perturbations from inflation; examples include modulated fluctuations [H, H| and the curvaton mechanism 
Q. These various scenarios all lead to similar predictions for the power spectrum. On the other hand, nongaussian 
statistics (such as the bispectrum) provide a powerful tool to observationally discriminate between different mecha- 
nisms for generating the curvature perturbation. In this review, which is based on , we will present a qualitatively 
new mechanism for generating cosmological perturbations during inflation. We discuss in detail the predictions of 
this new scenario for both the spectrum and nongaussianity of the primordial curvature fluctuations, showing how 
this new mechanism may be observationally distinguished from previous approaches. 



A. Nongaussianity from Inflation 

The possibility to discriminate between various inflationary scenarios has lead to a recent surge of interest in 
computing and measurin g nong aussian statistics. Although single field, slow roll models are known to produce 
negligible nongaussianity [1CH12| . there are now a variety of scenarios available in the literature which may predict 
an observable signature. Departures from gaussianity are often parametrized in the following form 

(( x ) = { g ( x ) + p NL [C g ( x )-(C g ( x ))] (1) 

where £(x) is the primordial curvature perturbation, Cg( x ) is a gaussian random field and /nl characterizes the degree 
of nongaussianity. The ansatz (p} is known as the "local" form of nongaussianity. 
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Although the local ansatz ([T]) has received significant attention, it is certainly not the only well-motivated model 
for a nongaussian curvature perturbation. For example, the nongaussian part of £(x) need not be correlated with the 
gaussian part. Consider a primordial curvature perturbation of the form 



where Fnl is some nonlinear (not necessarily quadratic) function and Xg( x ) is a gaussian field which is uncorrelated 
with Cg(x). Both (Q]) and {2J are local in position space, however, these two types of nongaussianity will have very 
different observational implications. The uncorrelated ansatz (JSJ for the primordial curvature perturbation can arise, 
for example, in models with preheating into light fields [HI, [HI]- (See also [HI, [l6j| for more discussion of nongaussianity 
from preheating and [l7| for another model where nongaussianity is generated at the end of inflation.) 

A useful quantity to consider is the bispectrum, B(ki, fc 2 , fca), which is the 3-point correlation function of the Fourier 
transform of the primordial curvature perturbation 



where ki = |ki|. The delta function appearing in ([3J reflects translational invariance and ensures that B(ki) depends 
on three momenta ki which form a triangle: ki + k2 + k3 = 0. Rotational invariance implies that B(ki) is symmetric 
in its arguments. 

If we assume the ansatz ([T]) for the primordial curvature perturbation then B(ki) has a very particular dependence 
on momenta; it peaks in the squeezed limit where one of the wave-numbers is much smaller than the remaining two 
(e.g. - ki <C k 2 , £3). Such a bispectrum is referred to as having a squeezed shape. However, other shapes of bispectrum 
are worth considering. A bispectrum is referred to as "equilateral" if it peaks when ki = k 2 = k 3 and "flattened" if 
it peaks when one of the wave-numbers is half the size of the remaining two (e.g. - 2ki = k 2 = k^). 

Without assuming any specific form for the primordial curvature perturbation, such as ([T]) or ([2]), one may charac- 
terize an arbitrary bispectrum ([3]) by specifying its shape, running and size [l8j |. As discussed above, the shape refers 
to the configuration of triangle on which B(ki) is maximal (squeeze, equilateral or flattened). The running of the 
bispectrum refers to how the magnitude of B(ki) depends on the overall size of the triangle. For example, in the case 
of scale invariant fluctuations, the bispectrum must scale as B{\k\ , Xk 2 , Xk^) = \~ & B{k\,k 2 ,ky,). Finally, the overall 
size of the bispectrum is often quantified by evaluating the magnitude of B(ki) on some fixed equilateral triangle. 
However, the skewness of the probability density function (defined later) might provide a better measure of the size 
of nongaussianity. 

Different types of nongaussian signatures are correlated with properties of the underlying inflation model. Let us 
first consider some examples with small running: 1 

1. A large bispectrum of local shape, along with iso-curvature effects, is associated with models where multiple 
fields are light (or otherwise dynamically important) during inflation. Examples include the curvaton mechanism 
[25I ] or models with turning points along the inflationary trajectory [2£|. The observational bound on local type 
nongaussianity, coming from the WMAP7 @ data, is -10 < /]^£ al < 74 [27j] at 95% confidence level. When 
combined with Large Scale Structure (LSS) data the bound becomes somewhat stronger: — 1 < /]^£ al < 65 [HJ. 

2. A large local bispectrum without any iso-curvature fluctuations can only be produced by nonlocal inflation 
models [2j| . For any single- field inflation model described by a local low-energy effective field theory, the results 
of [30I ] imply that the ratio of the 3-point correlation function to the square of the 2-point function must be 
of order the spectral tilt, in the squeezed limit. Hence, it has been argued that a large squeezed bispectrum 
must be associated with the presence of multiple light degrees of freedom, and hence iso-curvature effects. 
However, in [29j] it was shown single field nonlocal inflation models can produce a large squeezed bispectrum in 
the regime where the underlying scale of nonlocality is much larger than the Hubble scale during inflation. Such 
constructions evade the no-go theorem of [30| precisely because they violate the usual ass ump tion of cluster 
decomposition. Moreover, models of this type are not subsumed by the general analysis of [3l| since nonlocal 
field theories with infinitely many derivatives cannot be obtained in the regime of low-energy effective field 
theory. It is nevertheless sensible to study such constructions since they may be derived from ultra-violet (UV) 
complete frameworks, such as string field theory or p-adic string theory. See (32[ for details concerning the 
underlying consistency of nonlocal field theories and see (33| for a succinct review of nonlocal cosmology. 



C(a:) = C g (x) + F NL [ Xg (x)] 



(2) 



<Ck!Ck a Ck3> = (2^) 3 <5(k! + k 3 + k s )B(ki) 



(3) 



By "small" here we refer to any model where the running of the bispectrum is proportional to slow-variation parameters or arises due 
to loop effects. This does not necessarily mean that such running cannot lead to interesting observational signatures, see [2ll424ll . 
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3. A large equilateral bispectrum is typically associated with a small sound speed for the inflaton perturbations 
[HI], such as in Dirac-Born- Infold (DBI) inflation models [34[ . However, such a signature may also be ob- 
tained in multi-field gelaton [35[ or trapped inflation [36| models. The observational bound on equilateral type 
nongaussianity is —125 < //vx' 1 < 435 at 95% confidence level (37j . 

4. A large flattened bispectrum is associated with non- vacuum initial conditions p^l l38l - l40| . 2 

If we relax the assumption that the bispectrum is close to scale invariant then a much richer variety of nongaussian 
signatures is possible. For example, in models with sharp steps in the inflaton potential [4lT l42j the bispectrum is 
large only for triangles with a particular characteristic size. We will refer to such a signature as a localized nongaussian 
feature. Localized nongaussianities are not well constrained by current observation, but may be observable in future 
missions. 

Given the significant role that nongaussianity may play in discriminating between different models of the early 
universe, it is of crucial importance to explore and classify all possible consistent signatures for the bispectrum and 
other nongaussian statistics. Indeed, in this review we will describe a new kind of signature - uncorrelated nongaussian 
features - which is predicted in a variety of simple and well motivated models of inflation, but which has nevertheless 
been overlooked in previous literature. 



B. Inflationary Particle Production 

Recently, a new mechanism for generating cosmological perturbations during inflation was proposed @. This new 
mechanism, dubbed infra-red (IR) cascading, is qualitatively different from previous proposals (such as the curvaton or 
modulated fluctuations) in that it does not rely on the quantum vacuum fluctuations of some light scalar fields during 
inflation. Rather, the scenario involves the production of massive iso-curvature particles during inflation. These 
subsequently rescatter off the slow-roll condensate to generate bremsstrahlung radiation of light inflaton fluctuations 
(which induce curvature perturbations and temperature anisotropies in the usual manner). IR cascading can also 
be distinguished from previous mechanisms from the observational perspective: this new mechanism leads to novel 
features in both the spectrum and bispectrum. 

In principle, IR cascading may occur in any model where non-inflaton (iso-curvature) particles are produced during 
inflation. Models of this type have attracted considerable interest recently; examples have been studied where particle 
production occurs via parametric resonance [f| 0, HH, l43T - |48j , as a result of a phase transition [l3|, |49| - |55| , or otherwise 
[561 ] . Recent interest in inflationary particle production has been stimulated by various considerations: 

1. Particle production arises naturally in a number of microscopically realistic models of inflation, including ex- 
amples from string theory (36} and supersymmetric (SUSY) field theory [57| . In particular, inflationary particle 
production is a generic feature of open string inflation models Q, such as brane/axion monodromy (58l - l60| |. 

2. The energetic cost of producing particles during inflation has a dissipative effect on the dynamics of the inflaton. 
Particle production may therefore slow the motion of the inflaton, even on a steep potential. This gives rise 
to a new inflationary mechanism, called trapped inflation [36l . l6ll [621. which may circumvent some of the fine 
tuning problems associated with standard slow-roll inflation. See [361 ] for an explicit string theory realization 
of trapped inflation and [62j for a generalization to higher dimensional moduli spaces and enhanced symmetry 
loci. The idea of using dissipative dynamics to slow the motion of the inflaton is qualitatively similar to warm 
inflation [63[ and also to the variant of natural inflation [64[ proposed recently by Anber & Sorbo 56]. 

3. Observable features in the primordial power spectrum, generated by particle production and IR cascadin g, off er a 
novel example of the non-decoupling of high scale physics in the Cosmic Microwave Background (CMB) |43l[65|. 
In the most interesting examples, the produced particles are extremely massive for (almost) the entire history 
of the universe, however, their effect cannot be integrated out due to the non-adiabatic time dependence of the 
iso-inflaton mode functions during particle production. In (43| particle production during large field inflation 
models was proposed as a possible probe of Planck-scale physics. 



2 To our knowledge there is no explicit computation of the observational bound on flattened nongaussianity. In [38| a template (the 
enfolded model) was proposed. The analysis of [37|] is sufficiently general to study this shape, however, they do not explicitly place 
bounds on fff£ but instead constrain an alternative shape (the orthogonal model) which is a superposition of flattened and equilateral 
shapes. 
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In this article we study in detail the impact of particle production and IR cascading on the observable primordial 
curvature perturbations. In order to illustrate the basic physics we focus on a very simple and general prototype 
model where the inflaton, cj), and iso-infiaton, x, fields interact via the coupling 

Ant = -y(0-0o)V (4) 

We expect, however, that our results will generalize in a straightforward way to more complicated models, such as 
higher spin iso-inflaton fields or gauged interactions, wherein the physics of particle production and rescattering is 
essentially the same. Our result may also have implications for inflationary phase transitions, because spinodal decom- 
position can be interpreted as a kind of particle production and similar bi-linear interactions will induce rescattering 
effects. 

Scalar field interactions of the type (j4j have also been studied recently in connection with non-equilibrium Quantum 
Field Theory (QFT) |6p| . in particular with applications to the theory of preheating after inflation [67T - l7ll ] and also 
moduli trapping [6lL |62[ at enhanced symmetry points. Although our focus is on particle production during inflation 
(as opposed to during preheating, after inflation) some of our results nevertheless have implications for preheating, 
moduli trapping and also non-equilibrium QFT more generally. For example, in [6| analytical and numerical studies of 
rescattering and IR cascading during inflation made it possible to observe, for the first time, the dynamical approach 
to the turbulent scaling regime that was discovered in |72l [73j . 

Particle production during inflation in the model ((4]) leads to observable features in the primordial power spectrum, 
P{k). A number of recent studies have claimed evidence for localized features in P(k) that are incompatible with the 
simplest power-law model P(k) ~ k Tls ~ 1 [lil [73 - l86| . Although these observed features may simply be statistical 
anomalies (see, for example, [87] ) there remains the tantalizing possibility that they represent some new physics 
beyond the simplest slow roll model. Upcoming polarization data may play an important role in distinguishing these 
possibilities [54j |. In the meantime, it is interesting to determine the extent to which such features may be explained 
by a simple and well motivated model such as (U]). Moreover, because ^ is a complete microscopic model (as opposed 
to a phenomenological modification of the power spectrum) it is possible to predict a host of correlated observables, 
such as features in the scalar bispectrum and tensor power spectrum. Hence, it should be possible to robustly rule out 
(or confirm!) the possibility that some massive iso-curvature particles were produced during the observable portion 
of inflation. 

If detected, features from particle production and IR cascading will provide a rare and powerful new window into 
the microphysics driving inflation. This scenario opens up the possibility of learning some details about how the 
inflaton couples to other particles in nature, as opposed to simply reconstructing the inflaton potential along the 
slow roll trajectory. Moreover, due to the non-decoupling discussed above, features from particle production and IR 
cascading may probe new (beyond the standard model) physics at extraordinarily high energy scales. 

The outline of this paper is as follows. In section [II] we provide a brief, qualitative overview of the dynamics of 
particle production and IR cascading in the model ©. In section HiTl we study in detail this same dynamics using 
fully nonlinear lattice field theory simulations. In section IIVI we provide an analytical theory of particle production 
and IR cascading in an expanding universe. A complimentary analytical analysis, using second order cosmological 
perturbation theory, is provided in section [V] In section I VII we consider the observational constraints on inflationary 
particle production using a variety of data sets. In section I VIII we provide several explicit microscopic realizations 
of our scenario and study the implications of our observational constraints on models of string theory inflation, in 
particular brane monodromy. In section I Villi we quantify and characterize the nongaussianity generated by particle 
production and IR cascading. Finally, in section ITXl we conclude and discuss possible future directions. 



II. OVERVIEW AND SUMMARY OF THE MECHANISM 



In this section we provide a brief overview of the dynamics of particle production and IR cascading in the model 
((4]) and also summarize the resulting observational signatures. In the remainder of this article we will flesh out the 
details of this mechanism with analytical and numerical calculations. 

We consider the following model 



S = 



d A :> 



-fJ 



Ml 
~f R 



l -{d4>) 2 - V(ct>) - l -{d x f - 9 -{4> - 0o) V 



(•5) 



where R is the Ricci curvature constructed from the metric g^, 
usual, we assume a flat FRW space-time with scale factor a(t) 

ds 2 = g^ydx 11 dx v — —dt 2 



is the inflaton field and x 1S the iso-inflaton. As 



a 2 (t)dx 2 



(6) 
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and employ the reduced Planck mass M p = 2.43 x 10 18 GeV. We leave the potential V{(f>) driving inflation unspecified 
except to assume that it is sufficiently flat in the usual sense; that is e <C 1, <C 1 where 

M 2 /V'\ 2 V" 



are the usual slow roll parameters. 

Note that one might wish to supplement ([5]) by its supersymmetric completion in order to protect the flatness of 
the inflaton potential from large radiative corrections coming from loops of the x held. We expect that our results 
will carry over in a straightforward way to SUSY models and also to more complicated scenarios such as higher spin 
iso-inflaton fields and (possibly) inflationary phase transitions. 

The coupling 4>o) 2 X 2 m © is introduced to ensure that the iso-inflaton field can become instantaneously 

massless at some point <f> = <f>a along the inflaton trajectory (which we assume occurs during the observable range of 
e-foldings of inflation). At this moment x particles will be produced by quantum effects. 

Let us first consider the homogeneous dynamics of the inflaton field, 4>(t). Near the point (j> = (fio we can generically 
expand 

<P(t) ^(f> +vt (8) 

where v = (f)(0) and we have arbitrarily set the origin of time so that t — corresponds to the moment when <j) = <j>Q, 
(We are, of course, assuming that 0(0) 7^ 0.) The interaction ((4]) induces an effective (time varying) mass for the x 
particles of the form 

ml = g\4> - <i> ) 2 = kit 2 (9) 

where we have defined the characteristic scale 

K = v/flH (10) 

It is straightforward to verify that the simple expression (j9)) will be a good approximation for (iJ|t|) _1 < 0(e,rj) 
which, in most models, will be true for the entire observable 60 e-foldings of inflation. 

Note that, without needing to specify the background inflationary potential V((f>), we can write the ratio k*/H as 



h = I 9 

H y 2^l /2 



(11) 



where V^ 2 = 5 x 10 5 is the usual amplitude of the vacuum fluctuations from inflation. In this work we assume 
k± > H which is easily satisfied for reasonable values of the coupling g 2 > 10~ 7 . In particular, for g 2 ~ 0.1 we have 
fc*/F~30. 

The scenario we have in mind is the following. Inflation starts at some field value <p > 4>q and the inflaton rolls 
toward the point 4> — <j)Q. Initially, the iso-inflaton field is extremely massive m x 3> H and hence it stays pinned 
in the vacuum, x = 0, and does not contribute to super-horizon curvature fluctuations. Eventually, at t = 0, the 
inflaton rolls through the point <f> = 4>o where m x = and x particles are produced. To describe this burst of particle 
production one must solve for the following equation for the x-particle mode functions in an expanding universe 



Xk + 3Hxk + 



k 2 

- + kit 2 



Xk = (12) 



Equations of this type are well-studied in the context of preheating after inflation |68j and moduli trapping [61] . The 
initial conditions for (1121) should be chosen to ensure that the q-number field x is in the adiabatic vacuum in the 
asymptotic past (see sections IlIII and HVl for more details). In the regime fc* > H particle production is fast compared 
to the expansion time and one can solve (IT2|) very accurately for the occupation number of the created x particles 

n k = e - fe2 / fc * (13) 

Very quickly after the moment t — 0, within a time At ~ k^ 1 <C H~ , these produced x particles become non- 



-3 



relativistic (m x > H) and their number density starts to dilute as a 

Following the initial burst of particle production there are two distinct physical effects which take place. First, the 
energetic cost of producing the gas of massive out-of-equilibrium x particles drains energy from the inflaton condensate, 
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forcing 4> to drop abruptly. This velocity dip is the result of the backreaction of the produced x fluctuations on 
homogeneous condensate <j)(t). The second physical effect is that the produced massive \ particles rescatter off the 
condensate via the diagram Fig. [1] and emit bremsstrahlung radiation of light inflaton fluctuations (particles). 

Backreaction and rescattering leave distinct imprints in the observable cosmological perturbations. Let us first 
discuss the impact of backreaction. In Fig. [5] we plot the velocity dip resulting from the backreaction of the produced 
X particle on the homogeneous inflaton condensate <fi(t). From this figure we see that the quantity cj)/(H<fi) becomes 
large in the dip. This violation of slow roll is a transient effect; at late times the produced x particles become extremely 
massive and their number density dilutes as a~ 3 . 

One can understand the temporary slowing-down of the inflaton from an analytical perspective. Backreaction is 
taken into account using the mean-field equation 

+ 3ff0 + ^+ 5 2 (0-0 o )(x 2 ) =0 (14) 

where the vacuum average is computed following [U [6|| 

(X 2 ) = ^f^r (15) 

In equation (|14p we have implicitly assumed that the usual Coleman- Weinberg corrections to the inflaton potential 
have already been absorbed into V(4>), hence the vacuum average (x 2 ) should include only the effects of non-adiabatic 
particle production. (Here n x = J m£s n k ~ kl is the total number density of produced x particles and the factor 
a~ 3 reflects the usual volume dilution of non-relativistic matter.) In Fig. [2] we have plotted the solution of (fT4|) along 
with the exact result obtained from lattice field theory simulations, illustrating the accuracy of this simple treatment. 

Using the mean field approach, one finds that the transient violation of slow roll leads to a "ringing pattern" 
(damped oscillations) in the power spectrum P<p{k) — -^\5(j)k\ 2 of inflaton fluctuations [46T |, This ringing pattern is 
localized around wave-numbers which left the horizon at the moment when particle production occurred. The effect 
is very much analogous to Fresnel diffraction at a sharp edge. 

The second physical effect, rescattering, was considered for the first time in the context of inflationary particle 
production in 6]. Fig. [T] illustrates the dominant process: bremsstrahlung emission of long- wavelength Scf> fluctuations 
from rescattering of the produced x particles off the condensate. The time scale for such processes is set by the 
microscopic scale, A:" 1 , and is thus very short compared to the expansion time, H . Moreover, the production of 
inflaton fluctuations 5<p deep in the infra-red (IR) is extremely energetically inexpensive, since the inflaton is very 
nearly massless. The combination of the short time scale for rescattering and the energetic cheapness of radiating 
IR 5<f> leads to a rapid build-up of power in long wavelength inflaton modes: IR cascading. This effect generates a 
bump-like feature in the power spectrum of inflaton fluctuations, very different from the ringing pattern associated 
with backreaction. The bump-like feature from rescattering dominates over the ringing pattern from backreaction for 
all values of parameters. 

In the model (0 was studied using lattice field theory simulations, without neglecting any physical processes 
(that is to say that full nonlinear structure of the theory, including backreaction and rescattering effects, was accounted 
for consistently). However, this same dynamics can be understood analytically by solving the equation for the inflaton 
fluctuations 5(f) in the approximation that all interactions are neglected, except for the diagram Fig.[T] The appropriate 
equation is 

6$ + 3H5j> - %54> + = -g 2 \<j>(t) - M X 2 (16) 
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FIG. 2: \<j>\/(M p m) plotted against mt for g 2 = 0.1 (where m = is the effective inflaton mass). Time t = corresponds 
to the moment when <j) — cj>o and ^-particles are produced copiously. The solid red line is the lattice field theory result taking 
into account the full dynamics of rescattering and IR cascading while the dashed blue line is the result of a mean field theory 
treatment which ignores rescattering [4f|. The dot-dashed black line is the inflationary trajectory in the absence of particle 
creation. 



See Q for a detailed analytical theory. The solution of (fT6|) may be split into two parts: the solution of the homogeneous 
equation and the particular solution which is due to the source term. Schematically we have 

5<j>{t, x) = ^yacft.x) + <5<ft resc (t,x) (17) 

homogeneous particular 

The former contribution is the homogeneous solution which behaves as <5</> vac ~ H/(2tt) on large scales and, phys- 
ically, corresponds to the usual scale invariant vacuum fluctuations from inflation. The particular solution, S(j> resc , 
corresponds physically to inflaton fluctuations which are generated by rescattering. The abrupt growth of x inho- 
mogeneities at t = sources the particular solution <50 r esc, leading to the production of inflaton fluctuations which 
subsequently cross the horizon and become frozen. 

As mentioned earlier, rescattering generates a bump-like contribution to the primordial power spectrum of the 
curvature perturbations. To good approximation this may be described by a simple semi-analytic fitting function 

k \ ™ s 1 , /7re\ 3 / 2 / k N 



p(fc) = MW + MtJ UJ e ~ n7W (18) 

where the first term corresponds to the usual vacuum fluctuations from inflation (with amplitude A s and spectral 
index n s ) while the second term corresponds to the bump- like feature from particle production and IR cascading. 
The amplitude of this feature (^4ir) depends on g 2 , while the location (^ir) depends on 0rj. 

In Q the simple fitting function (|18p was used to place observational constraints on inflationary particle production 
using a variety of cosmological data sets. Current data are consistent with rather large spectral distortions of the 
type (|T8l . Features as large as O(10%) of the usual scale-invariant fluctuations from inflation are allowed, in the 
case that km falls within the range of scales relevant for CMB experiments. Such a feature corresponds to a realistic 
coupling g 2 ~ 0.01. Even larger values of g 2 are allowed if the feature is localized on smaller scales. In Fig. [3] we have 
illustrated the primordial power spectrum in the model ([5]) for a representative choice of parameters. We also plot 
the CMB angular Temperature- Temperature (TT) power spectrum for the same parameters. 

In 88| Large Scale Structure forecast constraints were considered for the model (TT51) . It was shown that, for 
km 0.1 Mpc -1 , the constraint on Am/A s will be strengthened to the 0.5% level by Planck or 0.1% including also 
data from a Square Kilometer Array (SKA). With a Cosmic Inflation Probe (CIP) similar constraints could be achieved 
for fciR. as large as 1 Mpc -1 . 

The prototype model ([5]) may be realized microscopically in a variety of different particle physics frameworks. 
In particular, particle production is a rather generic feature of open string inflation models [Tj where the inflaton, 
</>, has a geometrical interpretation as the position of some mobile D-brane. In this context the iso-inflaton, x, 
corresponds to a low-lying open string excitation which is stretched between the mobile inflationary brane and any 
other (spectator) branes which inhabit the compactification volume. If the inflationary and spectator branes become 
coincident during inflation, then the symmetry of the system is enhanced [6l| and some low-lying stretched string 
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FIG. 3: The top panel shows a sample bump in the power spectrum with amplitude Am, = 2.5 x 10~ which corresponds to a 
coupling g 2 ~ 0.01. The feature is located at A:ir = 0.01 Mpc -1 . This example represents a distortion of O(10%) as compared 
to the usual vacuum fluctuations and is consistent with the data at 2a. The bottom panel shows the CMB angular TT power 
spectrum for this example, illustrating that the distortion shows up mostly in the first peak. 



states will become instantaneously massless, mimicking the interaction (j4]) (see also [36(). An explicit realization of 
this scenario is provided by brane/axion monodromy models (58l - l60j . Our observational constraints on inflationary 
particle production may be used to place bounds on parameters of the underlying string model Q. 

The bump-like feature in P(k), illustrated in Fig. [3j must be associated with a nongaussian feature in the bispec- 
trum 0, i|. Indeed, it is evident already from inspection of equation (|T6|) that the inflaton fluctuations generated 
by rescattering are significantly nongaussian; the particular solution of (|16[) is bi-linear in the gaussian field x- The 
nongaussian signature from IR cascading is rather novel. The nongaussian part of £ is uncorrelated with the gaussian 
part. Moreover, the bispectrum B{ki) is very far from scale invariant; it peaks strongly for triangles with a character- 
istic size ~ fciR, corresponding to the location of the bump in the power spectrum ()18[) . The shape of the bispectrum 
therefore depends sensitively on the size of the triangle and is not well-described by any of the templates that have 
been proposed in the literature to date. 

The magnitude of this new kind of nongaussianity may be quite large. To quantify the effect it is useful to 
introduce the probability density function (PDF), P(£), which is the probability that the curvature perturbation has 
a fluctuation of size £. If we define the central moments of the PDF as 

(O = J CP(()dC (19) 
then a useful measure of nongaussianity is the dimensionless skewness of the PDF, defined by 

(C 3 )c 

where the subscript c indicates that only the connected part of the correlator should be included. The skewness S3 
encodes information about the bispectrum B{ki) integrated over all size and shape configurations and thus provides 
a meaningful single number to compare the nongaussianity of inflation models which may have very different shapes 
or running [2lj |. 
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If we choose g 2 ~ 0.01 (which is compatible with observation for all values of <po) then the model ([5]) produces the 
same value of S3 as a local model ([TJ with Jml — —53. This large value suggests that nongaussianity from particle 
production during inflation may be observable in future missions. 

Depending on model parameters, the nongaussian features predicted by the model ([5]) may lead to a rich variety 
of observable consequences for the CMB or Large Scale Structure (LSS). The phenomenology of this model is quite 
different from other constructions that have been proposed to obtain large nongaussianity from inflation. However, the 
underlying microscopic description ([3]) is extremely simple and, indeed, rather generic from the low-energy perspective. 
Explicit realizations of the interaction (|4]) have been obtained from string theory and SUSY. Moreover, in order to 
obtain an observable signature it was not necessary to fine-tune the inflationary trajectory or appeal to re-summation 
of an infinite series of high-dimension operators. 

III. NUMERICAL STUDY OF RESCATTERING AND IR CASCADING 

A. HLattice Simulations 

In this section we study numerically the creation of 5<f> fluctuations by rescattering of the produced x particles 
off the condensate 4>(t) in the model ([5]). To this end, we have written a new lattice field theory code, HLattice 
[89j |. for simulating the interactions of scalar fields in a cosmological setting. HLattice can be used to simulate the 
dynamics of any number of interacting scalar fields with arbitrary scalar potential and metric on field space [90l | . We 
solve the Klein-Gordon equations for the scalar field dynamics in an expanding FRW space-time and also solve the 
Friedmann equation self-consistently for the scale factor, a(t). Since the production of long wavelength 5<f> modes is so 
energetically inexpensive, a major requirement for successfully capturing this effect is respecting energy conservation 
to very high accuracy. HLattice conserves energy with an accuracy of order ~ 10 -8 , as compared to 10~ 3 — 10 -5 , 
which has obtained using previous codes such as DEFROST [9lJ or LATTICEASY [92|. A minimum accuracy of 
order 10 -4 is required for the problem at hand. 

The box size of our 512 3 simulations corresponds to a co-moving scale which is initially |^ ~ 3 times the horizon 
size H^ 1 , while fc* = 60^g H. We run our simulations for roughly 3 e-foldings from the initial moment t — when 
the x particles are produced, although a single e-folding would have been sufficient to capture the effect. For the sake 
of illustration, we have chosen the standard chaotic inflation potential V — m 2 4> 2 /2 with m = 10 _6 v87rM p for our 
numerical analysis. However, our results to not depend sensitively on the choice of background inflation model. 3 We 
have considered both fa = 2^/SttM p and 4>o = 3.2^/8ttM p and also three different values of the coupling constant: 
g 2 = 0.01,0.1,1. As expected, the coupling g 2 determines the magnitude of the effect while 4>o simply shifts the 
location of the power spectrum feature. For for this particular inflationary potential, the choice 4>o — S-2^/&nM p 
corresponds to putting the feature on scale slightly smaller than todays horizon. On the other hand, 4>o = 2\/%ttM p 
corresponds to placing the feature on scales much smaller than those probed by the CMB (we considered this case in 
order to be able to directly contrast our results with (46|V 

In order to capture the quantum production of x particles using classical lattice simulations, we start our numerical 
evolution very shortly after particle production has occurred, when the Xk modes are nearly adiabatic, but before 
any significant inflaton fluctuations have been produced. In practice, this corresponds to initializing the simulation 
at a time initial = 0{k^ ). The initial conditions for the modes Xk(t) are given by the usual Bogoliubov computation 
[UHII]. These are chosen to reproduce the occupation number = e~ 7!k l k * , while ensuring that the source term for 
the 6<p fluctuations is turned on smoothly at the initial time. As long as the initial conditions are chosen appropriately, 
our results are not sensitive to the choice of initial- 

At the initial time, the occupation numbers in the inflaton and iso-inflaton fluctuations are small. However, very 
quickly the massive x particles are diluted away by the expansion of the universe and the occupation number of 
the produced IR 5<j) fluctuations grows large compared to unity. Thus, classical lattice field theory simulations are 
sufficient to capture the late-time dynamics. (In the next section we will provide a quantum mechanical treatment of 
the dynamics of particle production and IR cascading during inflation, which will serve as an a posteriori justification 
for our classical lattice calculation.) 



The model independence of our result arises simply because all the dynamics of rescattering and IR cascading occurs within a single 
e-folding from the moment when </> = <po- Over such a short time it will always be a good approximation to expand <fi(t) = cf>0 + vt. 
Hence the dependence on the background dynamics arises only through v = 0(0) which is determined by the Hubble scale and the 
observed amplitude of curvature perturbations. This claim of model independence is born out by explicit analytical calculations in the 
next section. 
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Our approach is very similar to the methodology that has been employed successfully in studies of preheating after 
inflation for many years [9l|, 92]. In that case the initial fluctuations of the fields are chosen to reproduce the exact 
behaviour of the quantum correlation functions. The occupation numbers of the fields are small at the initial time. 
However, these grow rapidly as a result of the the preheating instability and classical simulations are sufficient to 
capture the late-time dynamics. 

B. Numerical Results 

We have studied the fully nonlinear dynamics of x particle production and the subsequent interactions of the 
produced x with the inflaton field in the model ([5]), as described above. We are interested in the power spectrum of 
the inflaton fluctuations 

^(fc) = ^W*l a (21) 

This contains a contribution coming from the usual quantum vacuum fluctuations from inflation that is close to the 
usual power- law form k na ~ l on large scales. Such a contribution would be present even in the absence of particle 
production and is not particularly interesting for us. In order to isolate the effects of rescattering we have subtracted 
off this component in figures |4l [5] and [6] In all cases we have normalized to the amplitude of the usual vacuum 
fluctuations from inflation, H 2 /(2ir) 2 . 

Fig. @] shows time evolution of the power in the inflaton fluctuations generated by rescattering, for three different 
time steps early in the evolution. This figure illustrates how multiple rescatterings leads to a dynamical cascading of 
power into the IR. To illustrate the magnitude of this effect, the horizontal yellow line corresponds to the amplitude 
of the usual vacuum fluctuations from inflation. For g 2 >, 0.06, the fluctuations from rescattering come to dominate 
over the vacuum fluctuations within a single e-folding. In Fig. [5] we illustrate how the magnitude of the spectral 
distortion depends on the coupling, g 2 . (The apparent change in the location of the feature for different values of g 2 
arises because we are plotting the power spectrum as a function of ln(fc/fc*) and k+ depends on g.) 

At late times, the IR portion of the power spectrum illustrated in Fig. [5] will remain fixed since the modes associated 
with these scales have gone outside the horizon and become frozen. On the other hand, the UV portion of this curve 
corresponds to modes that are still inside the horizon, hence we expect 6<frk ~ a" 1 and the UV tail of the power 
spectrum should damp as a -2 , due to the Hubble expansion. We observe precisely this behaviour in our lattice field 
theory simulations and this is illustrated in Fig. |6l which displays the dynamics of IR cascading over a much longer 
time scale. 

Within a few e-foldings from the time of particle production, the entire bump-like feature from IR cascading 
becomes frozen outside the horizon. At this point the fluctuations have become classical, large-scale adiabatic density 
perturbations and are observable in the present epoch (presuming that <p = 4>q occurs during the observable range of 
e-foldings) . In Fig. [3] we have illustrated this bump-like feature in both the primordial power spectrum and angular 
TT spectrum, for a representative choice of parameters. 

C. Backreaction Effects 

As discussed previously, the production of \ fluctuations at t = back-reacts on the homogeneous <p(t) causing a 
transient violation of slow roll. We can study this backreaction numerically, by averaging the inhomogeneous field 
4>(t, x) over the simulation box. The result is plotted in Fig. [2] We have also plotted the analytical solution of the 
mean- field equation (I14[) . showing that this agrees with the exact numerical result. 

The dynamics illustrated in Fig. [2] is easy to understand physically. The production of x particles at t = drains 
kinetic energy from the condensate and hence <f) must decrease abruptly. However, within a few e-foldings of the 
moment t = 0, the produced iso-inflaton particles become extremely massive and are diluted by the expansion as a -3 . 
At late times the inflaton velocity </> must tend to the slow roll value. Notice that the velocity including backreaction 
effects is not changed significantly, as compared to the usual slow roll result. This illustrates the energetic cheapness 
of particle production and IR cascading in the model (JS|). 

The transient violation of slow-roll illustrated in Fig. [5] is expected to induce a ringing pattern in the vacuum 
fluctuations from inflation [46| . This effect is accounted for automatically in our HLattice simulations. However, we 
would like to disentangle the effect of backreaction on the cosmological fluctuations from the effect of rescattering. 
This will be useful in order to compare the relative importance of different physical processes, and also to guide 
our analytical efforts in the next section. To this end, we consider the evolution of the curvature perturbation on 
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FIG. 4: The power spectrum of inflaton modes induced by rescattering (normalized to the usual vacuum fluctuations) as a 
function of ln(fc/fc*), plotted for three representative time steps in the evolution, showing the cascading of power into the IR. For 
each time step we plot the analytical result (the solid line) and the data points obtained using lattice field theory simulations 
(diamonds). The time steps correspond to the following values of the scale factor: a = 1.03,1.04,2.20 (where a = 1 at the 
moment when cj> — (fro). By this time the amplitude of fluctuations is saturated due to the expansion of the universe. The 
vertical lines show the range of scales from our lattice simulation. 




FIG. 5: The dependence of the power spectrum P^ on the coupling g 2 . The three curves correspond to P^, for g 2 = 0.01, 0.1, 1, 
evaluated at a fixed value of the scale factor, a — 2.20. We see that even for small values of g 2 the inflaton modes induced by 
rescattering constitute a significant fraction of the usual vacuum fluctuations after only a single e-folding. 
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FIG. 6: The power spectrum of inflaton modes induced by rescattering. (normalized to the usual vacuum fluctuations) as a 
function of ln(fc/fc*), plotted for three representative time steps in the late-time evolution. This figure illustrates the final stages 
of IR cascading; we see the peak of the bump-like feature slide to k ~ e _3 fc*, at which point the associated mode functions 
8(f)k have crossed the horizon and become frozen. At later times in the evolution the peak of the feature and also the IR tail 
(~ k 3 ) remain fixed. Modes associated with the UV end of the spectrum are still inside the horizon and continue to evolve as 
8<j)k ~ a -1 , which explains the damping of the k > e -2 fc* part of the spectrum. For each time step we plot the analytical result 
(the solid line) and the data points obtained using lattice field theory simulations (diamonds). 



co-moving hyper-surfaces, 1Z. In linear theory the equation for the Fourier modes TZk is well known 

z' 

n'j. + 2-74 + k 2 n k = (22) 

Here the prime denotes derivatives with respect to conformal time r = J — and z = a4>/H. Equation (|22|) is only 
strictly valid in the absence of entropy perturbations. However, in our case the \ held is extremely massive m 2 ^> H 2 
for nearly the entire duration of inflation, hence one may expect that direct iso-curvature contributions to 1Z are 
small. We have solved equation (|2"2"]l numerically. In order to take backreaction effects into account we compute the 
dynamics of z{t) = a(t)<j)(t) / H (i) by averaging over our HLattice simulation box. Next, we solve equation ([22]) given 
this background evolution and compute the power spectrum 

P n {k) = ^l|^ fe | 2 (23) 

The result is very close to the usual power-law form fc™ s_1 , with small superposed oscillations resulting from the 
transient violation of slow roll, see Fig. [7] In order to make the ringing pattern more visible, we have subtracted off 
the usual (nearly) scale-invariant result which would be obtained in the absence of particle production. For comparison, 
we also plot the bump-like feature from rescattering and IR cascading. This latter contribution was obtained using 
the results for P<j,{k) from the previous subsection and the naive formula TZ ~ £S<p (so that P n ~ (2eM 2 )- 1 P 4> ). 

From Fig. [7] we see that IR cascading has a much more significant impact on the observable curvature fluctuations 
than does backreaction. Indeed, for g 2 =0.1 the transient violation of slow- roll yields an order 10~ 2 correction to the 
vacuum fluctuations, while the correction from IR cascading is of order 10 . This dominance is generic for all values 
of the coupling. Thus, in developing an analytical theory of particle production during inflation, it is a very good 
approximation to completely ignore backreaction effects. 
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FIG. 7: A comparison of curvature fluctuations from different physical effects. The dashed blue line is the usual (nearly) 
scale invariant vacuum fluctuations from inflation. The red solid line is the bump-like feature induced by rescattering and IR 
cascading. The dotted blue line is the ringing pattern resulting from the momentary slowing-down of the inflaton (computed 
using the mean field approach of [461] ). The vertical lines show aH at the beginning of particle production and after ~ 3 
e-foldings. This figure clearly illustrates the dominance of IR cascading over backreaction effects. For illustration we have 
taken g 2 = 0.1, but the dominance is generic for all values of the coupling. 



IV. ANALYTICAL FORMALISM 



In the last section, we have studied particle production, rescattering and IR cascading using nonlinear lattice field 
theory simulations. In this section we will develop a detailed analytical theory, in order to understand those results 
from a physical perspective. These results were first presented in |8|. We consider, again, the model ([5]). The equations 
of motion that we wish to solve are 

-□0 + y'(0)+ 3 2 (0-0ok 2 = (24) 

-□x + 3 2 (0-0o) 2 X = (25) 

where □ = c^V^y is the covariant d'Alembertian. It will be useful to work with conformal time r, related to 
cosmic time t via adr = dt. In terms of conformal time the metric takes the form 

ds 2 = -dt 2 + a 2 (t)dx- dx 

= a 2 (r) [~dr 2 + dx ■ dx] (26) 



We denote derivatives with respect to cosmic time as / = dtf and with respect to conformal time as /' = d T f. 
The Hubble parameter H = a/ a has conformal time analogue % — a' /a. For an inflationary (quasi-de Sitter) phase 
(H = const) one has 

1 1 * = -V- (27) 



Ht 1 — € T 1 

to leading order in the slow roll parameter e <C 1. 

As discussed in sectionflll the motion of the homogeneous inflaton <j>(t) leads to the production of a gas of x particles 
at the moment t = when cj) = (f> . The first step in our analytical computation is to describe this burst of particle 
production in an expanding universe. Following the initial burst, both backreaction and rescattering effects take place. 
Our formalism will focus on the latter effect, which is much more important, and we provide only a cursory treatment 
of backreaction. 



A. Particle Production in an Expanding Universe 



The first step in our scenario is the quantum mechanical production of x-particles due to the motion of <fi. To under- 
stand this effect we must solve the equation for the x fluctuations in the rolling inflaton background. Approximating 
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6q + vt equation (|25j) gives 



X + ZHx - — X + k 4 J 2 X = 



(28) 



where fc* = J g\v\. We remind the reader that fc* 3> H for reasonable values of the coupling, see equation (|11[) . 

The flat space analogue of equation (|28[) is very well understood from studies of broad band parametric resonance 
during preheating [68j and also moduli trapping at enhanced symmetry points [61]. One does not expect this treatment 
to differ significantly in our case since both the time scale for particle production At and the characteristic wavelength 
of the produced fluctuations A are small compared to the Hubble scale: At ~ A -1 ~ k^ 1 <C H~ l , Hence, we expect 
that the occupation number of produced x particles will not differ significantly from the flat-space result (|13[) . at least 
on scales k >, H . Furthermore, notice that the x field is extremely massive for most of inflation 



1/2 



k I 



H 2 H 2 



(29) 



Since fc* ^ H it follows that m 2 , ^> H 2 , except in a tiny interval H\At\ ~ (ff/fc*) 2 which amounts to roughly 10 3 
e-foldings for g 2 ~ 0.1. Therefore, we do not expect any significant fluctuations of x to be produced on super-horizon 
scales k < H. 

Let us now consider the solutions of equation (|28[) . We work with conformal time r and write the Fourier transform 
of the quantum field x as 



x(t, x) 



d 3 k ^(t) 
(2tt)3/2 a (r) ' 



(30) 



Note the explicit factor of a 1 in (T5U|) which is introduced to give a canonical kinetic term. The q-number valued 
Fourier transform (r) can be written as 



£k(0=akXfe(r)+al k Xfc(T) 
where the annihilation/creation operators satisfy the usual commutation relation 

o k , 4'] =<5 (3) (k-k') 
and the c-number valued mode functions Xk obey the following oscillator-like equation 

Xk(r)+u;l(T)x k (T)=Q 

The time-dependent frequency is 



(31) 



(32) 



(33) 



2/ \ / 2 2 2 / \ " 

w fc (r) = k +a mJr) - — 



where to 2 (r) = g {4> — 4>o) — k*t (t) is the time-dependent effective mass of the x particles and 

t( T ) = J_ l n [ Zi 
1 ' H \Ht 



(34) 



(35) 



is the usual cosmic time variable. We have arbitrarily set the origin of conformal time so that r = —1/H corresponds 
to the moment when <f> — 4>q . 

In the left panel of Fig. [8] we have plotted a representative solution of (|33|) in order to illustrate the qualitative 
behaviour of the modes Xk- In the right panel of Fig.[S]we plot the occupation number rife of particles with momentum 
k, defined as the energy of the mode ^Ix'fel 2 + | w felXfc| 2 divided by the energy w k of each particle. Explicitly, we define 



Uk_ 

2 



' |2 



Ml 

w2 



+ IXfc| 5 



(36) 
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FIG. 8: The left panel illustrates the time dependence of the solutions \k of equation (|33[) for a representative choice of 
parameters. The oscillatory behaviour at early times represents the adiabatic initial condition. At t = the effective frequency 
(|34[) varies non-adiabatically and the fluctuations get a "kick" . The right panel plots the occupation number (|36[) for the same 
mode. No particles are present in the adiabatic in-going regime. This figure illustrates how the violations of adiabaticity at 
t = lead to production of \ particles. 



where the term — i comes from extracting the zero-point energy of the linear harmonic oscillator (see (68| for a review). 
From the left panel of Fig. [5]we see that, near the massless point t — 0, the fluctuations Xk get a "kick" and from the 
right panel we see that the occupation number n k jumps abruptly at this same moment. 

Let us now try to understand analytically the behaviour of the solutions of ([551) . At early times t <C — k" 1 , the 
frequency u> k varies adiabatically 



< 1 



(37) 



In this in-going adiabatic regime the modes Xk are not excited and the solution of (|33[) are well described by the 
adiabatic solution Xfe( r ) — /fc( T ) where 



fk(r) 



1 



\/2Wfe(7 



: exp 



i / dr'uj k (T) 



(38) 



We have normalized (138)) to be pure positive frequency so that the state of the iso-inflaton field at early times 
corresponds to the adiabatic vacuum with no x particles. (Inserting (|38[) into (|36[) one finds n k — for the adiabatic 
solution, as expected.) 

The adiabatic solution (1381) ceases to be a good approximation very close to the moment when <p = cf>o, that is at 
times \t\ < ■ In this regime the adiabaticity condition (|37[) is violated for modes with wave- number H <, k < fc* 
and x particles within this momentum band are produced. During the non-adiabatic regime we can still represent 
the solutions of (|33|) in terms of the functions f k (r) as 



Xk(r)=a k (T)f k {T)+h(T)r k {r) 



(39) 



This expression affords a solution of (|3"3")) provided the time-dependent Bogoliubov coefficients obey the following set 
of coupled equations 



P'kir) = 



2w fe (r) 
2w fc (r) 



exp 



exp 



-2i 



-2i 



drW(r') 
dr'uj k {T') 



Pk(r) 
Uk{r) 



(40) 
(41) 



The Bogoliubov coefficients are normalized as \u k \ 2 — \ j3 k \ 2 = 1 and the assumption that no x particles are present in 
the asymptotic past 4 fixes the initial conditions a k = 1, (3 k = for t — > —oo. This is known as the adiabatic initial 
condition. 



4 This assumption is justified since any initial excitation of \ would have been damped out exponentially fast by the expansion of the 
universe. 
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Prom the structure of equations (I40l4ip it is clear that violations of the condition (|37|) near t = leads to a rapid 
growth in the \/3 k \ coefficient. The time variation of /3 k can be interpreted as a corresponding growth in the occupation 
number of the x particles 

n k = \h\ 2 (42) 

At late times (t > k~ x ) adiabaticity is restored and the growth of n k = \(3 k \ 2 must saturate. By inspection of 
equations (|40I41|) we can see that the Bogoliubov coefficients must tend to constant values in the out-going adiabatic 
regime. Therefore, within less than an e-folding from the moment of particle production the solution Xk of equation 
P3p can be represented as a simple superposition of positive frequency f k modes and negative frequency /j* modes. 
Our goal now is to derive an analytical expression for the modes Xk which is valid in this out-going adiabatic region. 

Let us first study the adiabatic solution /^(t). If we focus on the interesting region of phase space, H < k < k+, 
then the adiabatic solution (f38|) is very well approximated by 

Mr) = 1 - e -^V) (43) 

where tir) is defined by (|35|) . It is interesting to note that equation (|43|) is identical to the analogous flat-space 
result Q, except for the factor of a -1 / 2 . Taking into account also the explicit factor of a -1 in our definition of the 
Fourier transform ([50)) we recover the expected large-scale behaviour for a massive field in de Sitter space, that is 
X ~ a -3 / 2 . This dependence on the scale factor is easy to understand physically, it simply reflects the volume dilution 
of non-relativistic particles: p x ~ m 2 % 2 ~ a~ 3 . 

Next, we seek an expression for the Bogoliubov coefficients ak, f3 k m the out-going adiabatic regime t £ k^ 1 . 
From (I40l4ip it is clear that the value of the Bogoliubov coefficients at late times can depends only on dynamics 
during the interval \t\ <, k^ 1 where the adiabaticity condition (|3"7| is violated. This interval is tiny compared to 
the expansion time and we are justified in treating a(r)&s a constant during this phase. Hence, it follows that the 
flat space computation of the Bogoliubov coefficients (6J, [68[ must apply, at least for scales k H. To a very good 
approximation we therefore have the well-known result 



a k = Vl + eF^V^ (44) 
(3 k ~ _i e -»*V(2^) (45) 

in the out-going adiabatic regime. Equation (|45[) gives the usual expression (|13[) for the co-moving occupation number 
of particles produced by a singe burst of broad-band parametric resonance: 

n k = \0 k \ 2 =e-"#rt (46) 

Finally, we arrive at an expression for the out-going adiabatic x modes which is accurate for interesting scales 
k± < k < H. Putting together the results (|43|) and (f39| along with the well-known expressions (|44I45[) we arrive at 

valid for t £ k* 1 . Equation (|47|) is the main result of this subsection. We will now justify that this expression is quite 
sufficient for our purposes. 

For modes deep in the UV, k £ fc*, our expression (|47[) . is not accurate. 5 However, such high momentum particles 
are not produced, the condition (|37|) is always satisfied for fc ^> fc*. Note that the absence of particle production deep 
in the UV is built into our expression (|47[) : as k — > oo this function tends to the vacuum solution Xk fk- 

Our expression (|47[) is also not valid deep in the IR, for modes k < H. To justify this neglect requires somewhat 
more care. Notice that, even very far from the massless point, t — 0, long wavelength modes k <C H should not be 
thought of as particle-like. The large-scale mode functions are not oscillatory but rather damp exponentially fast as 
X ~ a~ 3 / 2 . Hence, even if we started with some super- horizon fluctuations of x a t the beginning of inflation, these 



5 The expression H43I I for the adiabatic modes is not valid at high momenta where u)^ 
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would be suppressed by an exponentially small factor before the time when particle production occurs. Any super- 
horizon fluctuation generated near t = would need to be exponentially huge to overcome this damping. However, 
resonant particle production during inflation does not lead to exponential growth of mode functions. 6 

To verify explicitly that there is no significant effect for super-horizon fluctuations let us consider solving equation 
([28]) neglecting gradient terms. The equation we wish to solve, then, is 



dUa 3/2 x) + 



k\t 2 - -H 2 
4 



(a 3/2 x) = (48) 



(For simplicity we take e = for this paragraph, however, this has no effect on our results.) The solution of this 
equation may be written in terms of parabolic cylinder functions D u (z) as 

X(*,x)~ -I- (c x D 1 + gH 2 [{l+i)kJ] + C 2 D_ 1 9H2 [(-l + i)M]1 (49) 
a J/2 \ 3 H 7 ?" 3 1*5"' / 

For our purposes the precise values of the coefficients C\, C2 are not important. Rather, it suffices to note that for 
fc*|t| > 1 the function l|49p behaves as 

X(i,x) - |t|-V2 e -3«*/2 x [oscillatory] (50) 

This explicit large-scale asymptotics confirms our previous claims that the super-horizon fluctuations of \ damp to 
zero exponentially fast, as a~ 3 / 2 ~ e -3Ht/2^ ^ g discussed previously, this damping is easy to understand in terms of 
the volume dilution of non-relativistic particles. We can also understand the power-law damping that appears in (|50[) 
from a physical perspective. The properly normalized modes behave as o?^ 2 x ~ ^ while on large scales we have 
u>k ~ \ m x \ ~ Kl- Hence, the late-time damping factor £ -1 / 2 which appears in (150)) reflects the fact that the x particles 
become ever more massive as cj> rolls away from the point <pQ. 

Finally, it is straightforward to see that the function (j49j) does not display any exponential growth near t = 0. 
Hence, we conclude that there is no significant generation of super-horizon \ fluctuations due to particle production. 7 

In this subsection we have seen that the quantum production of x particles in an expanding universe proceeds very 
much as it does in flat space. This is reasonable since particle production occurs on a time scale short compared to 
the expansion time and involves modes which are inside the horizon at the time of production. 



B. Inflaton Fluctuations 



In section HV Al we studied the quantum production of x particles which occurs when <fi rolls past the massless point 
4> = 4>q. Subsequently, there are two distinct physical processes which take place: backreaction and rescattering. As 
we have argued in section IIII1 the former effect has a negligible impact of the observable spectrum of cosmological 
perturbations. Hence, we will not study this effect in any detail (see [i^,|4^,|4^] for analytical calculations). Instead 
we provide a cursory treatment of backreaction in Appendix A, in order to clear up some common misconceptions. 

In this subsection we study the rescattering of produced x particle off the inflaton condensate. The dominant 
process to consider is the diagram illustrated in Fig. [TJ corresponding to bremsstrahlung emission of 8<p fluctuations 
(particles) in the background of the external field. (There is also a sub-dominant process of the type XX 8<fi8(j) 
which is phase space suppressed.) Taking into account only the rescattering diagram illustrated in Fig.[T]is equivalent 
to solving the following equation for the q-number inflaton fluctuation 

8^ = -g 2 [0(i) - fa] x 2 (51) 



V 2 

8 2 + 3Hd t - ^ + to 2 
ar 



° In this regard our scenario is very different from preheating at the end of inflation. In the latter case the inflaton passes many times 
through the massless point m x = and there are, correspondingly, many bursts of particle production. After many oscillations of the 
inflaton field, the \ particle occupation numbers build up to become exponentially large and, averaged over many oscillations of the 
background, the \ mode functions grow exponentially. However, in our case case there is only a single burst of particle production at 
t = 0. The resulting occupation number J 13H is always less than unity and the solutions of J33H never display exponential growth. 

7 This is strictly true only in the linearized theory. It is possible that \ particles are generated by nonlinear effects such as rescattering. 
However, even such second order \ fluctuations will be extremely massive compared to the Hubble scale and must therefore suffer 
exponential damping a -3 / 2 on large scales. 
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where we have introduced the notation m 2 = for the infiaton effective mass. (Note that we are not assuming a 
background potential of the form to 2 </> 2 /2, only that ^ in the vicinity of the point <fi = (fig.) 

Equation (f5"Tj) may be derived by noting that (O gives an interaction of the form g 2 (<fi — 4>o)5(j)x 2 between the 
infiaton and iso-inflaton, in the background of the external field <fi(t). Equivalently, one may construct this equation 
by a straightforward iterative solution of (|24l) . 

We work in conformal time and define the q- number Fourier transform £k( r ) of the infiaton fluctuation analogously 
to (1301): 



(T ' xH /(2$7 



(2tt) 3 /2 a(r) 



(52) 



(To avoid potential confusion we again draw the attention of the reader to the explicit factor a 1 in our convention 
for the Fourier transform.) The equation of motion (|51l) now takes the form 



dl + k 2 + a 2 m 2 - — 



^(r)=-gkla(r)t(r) 



fX iX ( \ 
( 27r )3/2*k'Sk-k'W 



(53) 



The solution of (I53[) consists of two parts: the solution of the homogeneous equation and the particular solution which 
is due to the source. The former corresponds, physically, to the usual vacuum fluctuations from inflation. On the 
other hand, the particular solution corresponds physically to the secondary infiaton modes which are generated by 
rescattering. 



C. Homogeneous Solution and Green Function 



We consider first the homogeneous solution of f)53[) . Since the homogeneous solution is a gaussian field, we may ex- 
pand the q- number Fourier transform in terms of annihilation/creation operators 6k, 6^ and c-number mode functions 
fe (r) as 



tf(T) = ftkto(T) + &L k $S(r) 
Here the infiaton annihilation/creation operators 6k, 6^ obey 

bk,&£,] =<5 (3) (k-k') 

and commute with the annihilation/creation operators of the x~ field: 



a k ,6 



= 



(54) 



(55) 



(56) 



Using (|2"T)) and it is straightforward to see that the homogeneous infiaton mode functions obey the following 
equation 



d 2 T 4> k 



k 2 -\U 2 - 1 - 



4>k = o 



where we have defined 



u=-- V +e 



(57) 



(58) 



The properly normalized mode function solutions are well known and may be written in terms of the Hankel function 
of the first kind as 



</>k(r) 



-rHi l \-kr) 



(59) 



This solution corresponds to the usual quantum vacuum fluctuations of the infiaton field during inflation. 

In passing, let us compute the power spectrum of the quantum vacuum fluctuations from inflation. Using the 
solutions (l59l) we have 
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on large scales k <C aH. The explicit factor of a 2 in (|60|) appears to cancel the a 1 in our definition of the Fourier 
transform (|52[) . The spectral index is 



1 = 3 - 2v = 2r/ - 2e 



(61) 



using ((55)1 . 

Given the solution (|59l) of the homogeneous equation, it is now trivial to construct the retarded Green function for 
equation (|53j) . This may be written in terms of the free theory mode functions (|59|) as 



G k (r-T') = ie(r-r')[ ^(t)^(t') - ^(r)^(r') ] 



ITT 



(62) 



D. Particular Solution: Rescattering Effects 

We now consider the particular solution of (|53|) . This is readily constructed using the Green function (|62|) as 

&M=--(0L j dT'd*k'G k {T-Tl)a{T')t{T')ZUl-*{T') 

Notice that this particular solution is statistically independent of the homogeneous solution (|54|) . In other words, the 
particular solution can be expanded in terms of the annihilation/creation operators a^, associated with the \ field, 

whereas the homogeneous solution is written in terms of the annihilation/creation operators foki^k associated with 
the inflaton vacuum fluctuations. These two sets of operators commute with one another. 

We will ultimately be interested in computing the n-point correlation functions of the particular solution (|63|) . For 
example, carefully carrying out the Wick contractions, the connected contribution to the 2-point function is 

(^^M)-^ (3) (ki + k 2 ) 

X J dT'dT"a{T')a{T")t{T')t{T")G kl {T-T')G k2 {T-T") 

x J d 3 fc' Xfcl -^(r')Xfe 1 -^(r")x fc '(r')x^(r") (63) 

The power spectrum of 8(f) fluctuations generated by rescattering is then defined in terms of the 2-point function in 
the usual manner 

(ZimtiT)) = ^ (3) (k + k')^a(r) 2 P— (64) 

(The explicit factor of a 2 in the definition (|64[) appears to cancel the factor of a -1 in our convention for Fourier 
transforms (f52| .) 

The total power spectrum is simply the sum of the contribution from the vacuum fluctuations (|60|) and the contri- 
bution from rescattering ((M)) : 

p <j> (k) = p;^(k) + p^ c (k) (65) 

There are no cross-terms, owing to the fact and commute. 



E. Renormalization 



We now wish to evaluate the 2-point correlator ((63)) . In principle, this is straightforward: first substitute the result 
(j4"T)l for the \k modes and the result ((6"2")) for the Green function into (|63p . next evaluate the integrals. However, 
there is a subtlety. The resulting power spectrum is formally infinite. Moreover, the 2-point correlation function 
(|63p receives contributions from two distinct effects. There is a contribution from particle production, which we 
are interested in. However, there is also a contribution coming from quantum vacuum fluctuations of the \ field 
interacting non-linearly with the inflaton. The latter contribution would be present even in the absence of particle 
production, when a k = 1, /3fc = 0. 
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In order to isolate the effects of particle production on the inflaton fluctuations, we would like to subtract off the 
contribution to the 2-point correlation function (|63p which is coming from the quantum vacuum fluctuations of x- This 
subtraction also has the effect of rendering the power spectrum (|64|) finite, since it extracts the usual UV divergent 
contribution associated with the Minkowski-space vacuum fluctuations. 

As a step towards renormalizing the 2-point correlation function of inflaton fluctuations from rescattering (|63p , let 
us first consider the simpler problem of renormalizing the 2-point function of the gaussian field \- We defined the 
renormalized 2-point function in momentum space as follows: 

(&(tl)&(k)>ren - (&(*l)&(*2)) ~ 

In (f66|) the quantity (ti)£% (ia))in is the contribution which would be present even in the absence of particle 
production, computed by simply taking the solution (|39[) with a k = 1, = 0. Explicitly, we have 

<&(*l)&(*2)>in = S^(k 1 + k 2 )f kl ( tl )fl(t 2 ) (66) 

where f k are the adiabatic solutions (|38l) . 



To see the impact of this subtraction, let us consider the renormalized variance for the iso-inflaton field, (x 2 ) 
Employing the prescription (|66[) we have 



(X 2 (r,x)) r 



d 3 k 



(27r) 3 a 2 (r) 



(x 2 (t,x))-<5 m (67) 



n x = I 4^ n k ~ k l (69) 



where 5m is the contribution from the Coleman- Weinberg potential. This proves that our prescription reproduces the 
scheme advocated in [6l| . The renormalized variance (|S7| is finite and may be computed explicitly using our solutions 
07]). We find 

<X 2 (;,x)) rcn -^-!- (68) 

g\(p-<po\ 

where 

„. = r 

(2tt) 3 

is the total co-moving number density of produced x particles. The result (ff55|) was employed in [46[ to quantify the 
effect of backreaction on the inflaton condensate in the mean field treatment ([Tff]) . Hence, the renormalization scheme 
(|66| was implicit in that calculation also. 

At the level of the 2-point function, our renormalization scheme is tantamount to assuming that Coleman- Weinberg 
corrections are already absorbed into the definition of the inflaton potential, V((j>). In general, such corrections 
might steepen V{<f>) and spoil slow-roll inflation. Here, we assume that this problem has already been dealt with, 
cither by fine-tuning the bare inflaton potential or else by including extended SUSY (which can minimize dangerous 
corrections). See also [6l[ for a related discussion. Note also that our renormalization procedure is equivalent to the 
quasi-particle normal ordering scheme described in [93[ . 

Having established a scheme for remormalizing the 2-point function of the gaussian field Xj it is now straightforward 
to consider higher order correlation functions. We simply re- write the 4-point function as a product of 2-point functions 
using Wick's theorem. Next, each Wick contraction is renormalized as (|66p . Applying this prescription to (1631) amounts 
to 

«t(r)£(T)U = g^(k 1+ k 2 ) 

x J dr'dT"t(T')t(T")a(T')a(r")G kl (r - r')G k2 (r - r") 

x ( d 3 k> [x kl - k '(r')xt- k/ (r") - f kl -k'(r')fl-k'(r")] fok'COxMO " M*UMO] (70) 



where fk(i) are the adiabatic solutions defined in ([55)) . 
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F. Power Spectrum 



We are now in a position to compute the renormalized power spectrum of inflation fluctuations generated by 
rescattering, P£ esc (fc). We renormalize the 2-point correlator of the inflaton fluctuations generated by rescatter 
according to ([TO]) and extract the power spectrum by comparison to We have relegated the technical details to 

appendix B and here we simply state the final result 

„2pjL r -7rfc 2 /(2fcJ) 



P^ sc (k) 



16tt 5 

-7Tk 2 /(4k'i 



2V2 
1 

2V2 C 



-37r/c J /(8fc 2 ) 



8 'V / 2 ie -7rfe 2 /(3^) , 4 ^ c -37rfc 2 /(5fc 



3a/3 



5\/5 



7 2 (fc,r) 2 -Re [Ji(fc,r)]) 
Im [h(k,T)I 2 (k,T)} 



(71) 



where the functions Ii , / 2 are the curved space generalization of the characteristic integrals defined in [6| . Explicitly 
we have 



h(k,T) 

h{k,r) 



4t / dr>G k {T~T>y k2 * t2 ^ 

a ( T ) J 

dT'G k (T~T') 



a(r) 



(72) 
(73) 



The characteristic integral I 2 can be evaluated analytically, however, the resulting expression is not particularly 
enlightening. Evaluation of the integral Zj requires numerical methods. More details in Appendix B. Equation (|71|) 
is the main result of this section. 



G. Comparison to Lattice Field Theory Simulations 

In section IIIII the results of our analytical formalism were plotted alongside the out- put of fully nonlinear HLattice 
simulations. It is evident from Figs. 21 [5] and [6] that the agreement between these approaches is extremely good, even 
very late into the evolution and in the regime g 2 ~ 1. The consistency of perturbative quantum field theory analytics 
and nonlinear classical lattice simulations provides a highly nontrivial check on our calculation. 



H. The Bispectrum 



So far, we have shown how to compute analytically the power spectrum generated by particle production, rescatter- 
ing and IR cascading during inflation in the model (O . We found that IR cascading leads to a bump- like contribution 
to the primordial power spectrum of the inflaton fluctuations. However, this same dynamics must also have a nontriv- 
ial impact on nongaussian statistics, such as the bispectrum. Indeed, it is already evident from our previous analysis 
that the inflaton fluctuations generated by rescattering may be significantly nongaussian. From the expression (|63[) 
we see that the particular solution (due to rescattering) is bi-linear is the gaussian field x- 

We define the bispectrum of the inflaton field fluctuations in terms of the three point correlation function as 

W) = (2^) 3 a 3 (r)<5(k 1 + k 2 + k 3 )B^) (74) 

The factor a 3 appears in (17411 to cancel the explicit factors of a -1 in our convention l]52[) for the Fourier transform. 
It is well-known that the nongaussianity associated with the usual quantum vacuum fluctuations of the inflaton is 
negligible [Iol - fl2T | , therefore, when evaluating the bispectrum (|74| we consider only the particular solution (|63|) which 
is due to rescattering. Carefully carrying out the Wick contractions, we find the following result for the renormalized 
3-point function 

J d 3 p [xki-pirijXki-pfo) - fk 1 -p(n)fk 1 -p('n2)\ [Xk 3 +p(T2)xt 3 +p( T 3) - fk 3 + P (T2)ft 3+p (T 3 )] [x p (ti)x*(t 3 ) - / P (ri)/* 
+ (k 2 O h) 
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where the modes \k are defined by (|39|) and fk are the adiabatic solutions (|38|) . On the last line of (j75[) we have 
labeled schematically terms which are identical to the preceding three lines, only with fc 2 and k^ interchanged. One 
may verify that this expression is symmetric under interchange of the momenta fej by changing dummy variables of 
integration. 

It is now straightforward (but tedious) to plug the expressions (|43"]) and (|47l) into (|75|) and evaluate the integrals. 
This computation is tractable analytically because the time and phase-space integrals decouple. The bispectrum is 
then extracted by comparison to (|74l) . This computation is carried out in detail in Q where we have shown that 
B,f>(ki) peaks only over triangles with a characteristic size, corresponding to the location of the bump in the power 
spectrum. This result is easy to understand on physical grounds, all the dynamics of rescattering and IR cascading 
take place over a very short time near the moment t = 0. Hence, the effect of this dynamics on the primordial 
fluctuations must be limited to scales leaving the horizon near the time when particle production occurs. 

We will provide a cursory discussion of the bispectrum B^^ki) in section [Villi when we discuss nongaussianity from 
particle production during inflation. 



I. Inclusion of a Bare Iso-Inflaton Mass 



In passing, it may be interesting to consider how the analysis of this section is modified in the case that our 
prototype action ([5]) is supplemented by a bare mass term for the iso-inflaton field of the form 5C = — \[x 2 x 2 ■ Thus, 
in place of ([5]) suppose we consider the model 

£ = ~\{d<t>) 2 - V&) - \{d X f - \f X 2 ~ y (0 - to) V (76) 

Now the x particles do not become massless at the point <p = 4>o, but rather the effective mass-squared 

m\ =^+.g 2 W-0o) 2 (77) 

reaches a minimum value /i 2 (which we assume to be positive) . Such a correction may arise due to a variety of effects 
and will reduce the impact of particle production and IR cascading on the observable cosmological fluctuations. 

In [8] we have carefully calculated the spectrum and bispectrum for the model (76) and shown explicitly how these 
differ from the analogous results in the case /i 2 = 0. Rather than repeat this analysis, let us instead briefly consider 
how the additional bare mass term in (|76[) alters the dynamics of particle production from a more heuristic perspective. 
The iso-inflaton fluctuations now obey the equation 

X + 3H X - ^X + [p 2 + Kt 2 ] X = (78) 

rather than (|28|) . This equation was solved in (olj in the regime where particle production is fast compared to the 
expansion time. 8 The occupation number of produced x particles is 

ii t = re- rf2/l * (79) 
which differs from our previous result (fT3|) by the suppression factor 

r = e -^ 2 /^ < i ( 80 ) 

Therefore, the effect of the inclusion of a bare mass for the iso-inflaton is to suppress the number density of produced 
X particles by an amount r. This suppression reflects the red uced pha se space of produced particles: the adiabaticity 
condition |w^,/a; 2 | <C 1 is violated only for modes with k < \Jk1 — /i 2 . 

The reduction of rifc translates into a suppression for the n-point correlation functions of the iso-inflaton. For 
example, the renormalized variance (x 2 )ren cx J d 3 kn,k is suppressed by a factor of r. The power spectrum of inflaton 
fluctuations generated by rescattering, Pij csc , is proportional to the 4-point correlator of x and hence picks up a 



In the opposite regime, which corresponds to a fine-tuned coupling g 2 <C 10 — 7 , the iso-inflaton will be light for a significant fraction 
of inflation. In that case the theory 10 must be considered as a multi-field inflation model and one can no longer consistently assume 
(x) = 0- In other words, relaxing the assumption of fast particle production significantly changes the scenario under consideration and 
we do not pursue this possibility any further. 
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suppression factor of order r 2 . Similarly, the bispectrum B<p is proportional to the 6-point correlator of x and must 
be reduced by a factor of order r 3 . The condition 

,j 2 « kl (81) 

is equivalent to r = 1 and ensures that the addition of a bare iso-inflaton mass will have a negligible impact on any 
observable. 

For models obtained from string theory or super-gravity (SUGRA), it is natural to expect fi of order the Hubble 
scale during inflation [95T - [97l | . In the context of SUGRA, the finite energy density driving inflation breaks SUSY 
and induces soft scalar potentials with curvature of order Vj," ft ~ /i 2 ~ H 2 [96]. In the case of string theory, many 
scalars are conformally coupled to gravity [97[ through an interaction of the form 5C = —j^R\ 2 where the Ricci 
scalar is R ~ H 2 during inflation. More generally, any non-minimal coupling SC = — |i?x 2 between gravity and the 
iso-inflaton will induce a contribution of order H to the effective mass of x, as long as £ = 0(1). In all models where 
/i 2 ~ H 2 the condition (fgTj) is satisfied for reasonable values of the coupling g 2 > 10~ 7 , see equation (fTTj) . Thus, we 
expect that corrections of the form SC = — \^ 2 X 2 w iU n °t alter our results in a wide variety of well-motivated models. 
This heuristic expectation was verified explicitly in Q . 



V. COSMOLOGICAL PERTURBATION THEORY 



In section IIVI we developed an analytical theory of particle production and IR cascading during inflation which 
is in very good agreement with nonlinear lattice field theory simulations. However, this formalism suffers from a 
neglect of metric perturbations and, consequently, we were unable to rigorously discuss the gauge invariant curvature 
perturbation £. 9 Hence, the reader may be concerned about gauge ambiguities in our results. In this section we 
address such concerns, showing that metric perturbations may be incorporated in a straightforward manner and that 
their consistent inclusion does not change our results in any significant way. We will do so by showing explicitly that, 
with appropriate choice of gauge, equations (|51l) and (|2"5|) for the fluctuations of the inflaton and iso-inflaton still hold, 
to first approximation. We will also go beyond our previous analysis by explicitly showing that in this same gauge 
the spectrum of the curvature fluctuations, Pq, is trivially related to the spectrum of inflaton fluctuations, P^. 

To render the analysis tractable we would like to take full advantage of the results derived in the last section. To 
do so, we employ the Seery et al. for malis m for working directly with the field equations (98j and make considerable 
use of results derived by Malik in [99|, llOOl ] . (Note that our notations differ somewhat from those employed by Malik. 
The reader is therefore urged to take care in comparing our formulae.) 

We expand the inflaton and iso-inflaton fields up to second order in perturbation theory as 

0(r,x) = 0(T) + M(r,x) + ^ 2 0(r,x) (82) 

X(r,x) = 5ix(t,x) + -<J 2 x(r,x) (83) 

The perturbations are defined to average to zero (5 n <t>) = (S n x) — so that (<j>(t, x)) = <j>(t) and (x(i, x)) = 0. (The 
condition (x) = is ensured by the fact that m x 3> H for nearly the entire duration of inflation.) 

We employ the flat slicing and threading throughout this section. With this gauge choice the perturbed metric 
takes the form 

5 oo = -a 2 (l + 2^+^ 2 ) (84) 



9oi = a d, 



B, + -B 2 



2 

2 



(85) 



gij = a Sij (86) 

so that spatial hyper-surfac es ar e flat. Note also that in this gauge the field perturbations 6 n 4>, 8 n \ coincide with the 
Sasaki-Mukhanov variables |l0l| at both first and second order. 

This perturbative approach, of course, neglects the momentary slow-down of the inflaton background due to back- 
reaction. However, we have already shown in subsection IIII Gl that backreaction has a tiny impact on the observable 
cosmological perturbations (see also [6J). 



This variable is related to the quantity 1Z defined in subsection IIII Cl as £ = —72- on large scales. 
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A. Gaussian Perturbations 

In [9^| Malik has derived closed-form evolution equations for the field perturbations 5 n (f>, S n X a t both first (n = 1) 
and second [n = 2) order in perturbation theory. Let us first study the gaussian perturbations. The closed-form 
Klein-Gordon equation for Si4> derived in [99] can be written as 



<Si0" + ZHSkP - V 2 <5i^ 



2 2 
a m 



8x<j) = 



(87) 



Following our previous analysis we expand the first-order perturbation in terms of annihilation/creation operators as 



<Sl<K*)X) 



d 3 k 
(2tt) 3 / 2 



: a(r) ' 



.ik'X 



h.c. 



where h.c. denotes the Hcrmitian conjugate of the preceding term and we draw the attention of the reader to the the 
explicit factor of a -1 in our definition of the Fourier transform. Working to leading order in slow roll parameters we 
have 



1 



k z + — (-2 + 37? - 9e) 

T 



(89) 



This equation coincides exactly with (|57[) and the properly normalized solutions again take the form (|59p . The only 
difference is that the order of the Hankcl function, v, is now given by 



(90) 



rather than by equation (|58j) . The power spectrum of the gaussian fluctuations is, again, given by (|60)) . The correction 
to the order of the Hankcl function v translates into a correction to the spectral index: instead of (|61|) we now have 



n s — 1 = 2r] — 6e 



(91) 



which is precisely the standard result [94j . 

Thus, as far as the quantum vacuum fluctuations of the inflaton are concerned, the only impact of consistently 
including metric perturbations is an 0(e) correction to the spectral index n s . 

Let us now turn our attention to the first order fluctuations of the iso-inflaton. The closed-form Klein-Gordon 
equation for 8\x derived in [99j can be written as 



£ix" + mSix' - V 2 5 lX + a 2 fe^ 2 (r)<5i0 = 



(92) 



This coincides exactly with equation (f28|) . which we have already solved. The fact that linear perturbations of \ do 
not couple to the metric fluctuations follows from the condition (x) = 0. 



B. Nongaussian Perturbations 

Now let us consider now the second order perturbation equations. The closed-form Klein-Gordon equation for <52<^ 
derived in [99( can be written as 

8 2 <j)" + 2U8 2 <j)' - ^ 2 5 2 <j) + 

As usual, the left-hand-side is identical to the first order equation ([57]) while the source term J is constructed from 
a bi-linear combination of the first order quantities 8\(j) and S\\. In order to solve equation (|93[) we require explicit 
expressions for the Green function Gk and the source term J. The Green function is trivial for the case at hand; it is 
still given by our previous result (|5"5]) . provided one takes into account the fact that the order of the Hankel functions 
v is now given by (1901) . rather than (|58[) . In other words, the Green function for the nongaussian perturbations (|93[) 
differs from the result obtained neglecting metric perturbations only by 0(e) corrections. 



2_2 

a m 



M, : 



5 2 cf> = J(t, x) 



(93) 
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Next, we would like to consider the source term, J, appearing in (j93J). Schematically, we can split the source into 
contributions bi-linear in the gaussian inflaton fluctuation Sid and contributions bi-linear in the iso-inflaton Six- 

J=J (t> + J x (94) 

The contribution would be present even in the absence of the iso-inflaton. These correspond, physically, to the 
usual nongaussian corrections to the inflaton vacuum fluctuations coming from self-interactions. This contribution 
to the source is well-studied in the literature and is known to contribute negligibly to the bispectrum [98]. Thus, in 
what follows, we will ignore J^. 

On the other hand, the contribution J x appearing in (|94p depends only on the iso-inflaton fluctuations 5\X- This 
contribution can be understood, physically, as generating nongaussian inflaton fluctuations 52<fi by rescattering of the 
produced \ particles off the condensate. Hence, the contribution J x may source large nongaussianity and is most 
interesting for us. It is straightforward to compute J x explicitly for our model using the general results of [99(. We 
find 

J x = -2a 2 g 2 (6-6 )(5 lX ) 2 

-a 2 g 2 {6 - M 2 (Si X ) 2 - i(v^x) 2 - \(5i X ') 2 

+V- 2 (d,(<5 lX )V 2 ,9 l (<Six) + V 2 (<5 lX )V 2 (<5 lX ) + 5 lX "V 2 5 lX + (VS lX ') 2 ) ] (95) 

where the upper sign is for 6' > 0, the lower sign is for 6' < 0. Notice that the contributions to J x on the third and 
fourth line of (|95j) contain the inverse spatial Laplacian V -2 and are thus nonlocal. Th ese t erms all contain at least 



V2e 
M p 



as many gradients as inverse gradients and hence the large scale limit is well-defined. In [102j is was argued that these 
terms nearly always contribute negligibly to the curvature perturbation on large scales. 

Let us now examine the structure of the iso-inflaton source J x , equation (|9"5|) . The first line of ([§5"]) goes like 
a 2 g 2 (d — (po)(5ix) 2 ■ This coincides exactly with the source term in equation (|51[) which was already studied in section 
IIVI On the other hand, the terms on the second, third and fourth lines of ([95]) are new. These represent corrections to 
IR cascading which arise due to the consistent inclusion of metric perturbations. We will now argue that these "extra" 
terms are negligible as compared to the first line. If we denote the energy density in gaussian iso-inflaton fluctuations 
as p x ~ to 2 (Six) 2 then, by inspection, we see that the first line of ([95]) is parametrically of order p x /\4> — <fco\ while 
the remaining terms are or order yfep x /M p . Hence, we expect the first term to dominate for the field values 6 = 4>q 
which are relevant for IR cascading. This suggests that the dominant contribution to J x is the term which we have 
already taken into account in section IIVI 

Let us now make this argument more quantitative. Inspection reveals that the only "new" contribution to (|95[) 
which has any chance of competing with the "old" term a 2 g 2 (6 — (f>o)(5ix) 2 is the one proportional to y/ea 2 g 2 (6 — 
6q) 2 {8ix) 2 1 M p (the first term on the second line). This new correction has the possibility of becoming significant 
because it grows after particle production, as 6 rolls away from 6$. This growth, which reflects the fact that the 
energy density in the x particles increases as they become more massive, cannot persist indefinitely. Within a few 
e-foldings of particle production the iso-inflaton source term must behave as J x ~ a~ 3 , corresponding to the volume 
dilution of non-relativistic particles. Hence, in order to justify the analysis of section HVl we must check that the term 

■/new- j^a 2 g 2 (6~do) 2 (6 lX ) 2 (96) 

does not dominate over the term which we have already considered 

Joid- aV(0-0o)(<$ix) 2 (97) 
during the relevant time HAt = 0(1) after particle production. It is straightforward to show that 

J oW M p 1 M p H 1 1 1 



-'"now %/e <t> - 4>o 4>y/e N e N 



(98) 



where N = Ht is the number of e-foldings elapsed from particle production to the time when IR cascading has 
completed. Hence, N = 0(1) and we conclude that the second, third and fourth lines of (l95l are (at least) slow roll 
suppressed as compared to the first line. 

In summary, we have shown that consistent inclusion of metric perturbations yields corrections to the inflaton 
fluctuations 66 which fall into two classes: 
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1. Slow-roll suppressed corrections to the inflaton vacuum fluctuations 5\4> (these amount to changing the definition 
of v in the solution ([53])). These corrections have two physical effects. First, they yield an 0(e) correction to 
the spectral index. Second, they modify the propagator Gk by an 0(e) correction. 

2. Corrections to the source J for the nongaussian inflaton perturbation 52<fi. These corrections are the second, 
third and fourth lines of (195[1 ) which, as we have seen, are slow roll suppressed. 

It should be clear that neither of these corrections alters our previous analysis in any significant way. 

C. Correlators 

So far, we have shown that a consistent inclusion of metric perturbations does not significantly alter our previous 
results for the field perturbations. Specifically, Six is identical to our previous solution of equation (|28[) for the iso- 
inflaton, while 5\4> coincides with the homogeneous solution of equation (I51[) , up to slow-roll corrections. At second 
order in perturbation theory, we have seen that 

<5 2 = J d A x'G{x - x')J x {x') + [(<5i</>) 2 ] 

To leading order in slow roll J x = ~2a 2 g 2 (<fi — 4>o){$ix) 2 an d the first term coincides with our previous result for 
the particular solution of equation (|51| . The terms of order (<5i0) 2 represent nongaussian corrections to the vacuum 
fluctuations from inflation (coming from self-interactions of 6(j> and the nonlinearity of gravity). These would be 
present even in the absence of particle production, and are known to have a negligible impact on the spectrum and 
bispectrum [98| . 

We are ultimately interested in the connected n-point correlation functions of 6<j). For example, the 2-point function 
((6(f))' 2 ) get a contribution of the form ((di<p) 2 } which gives the usual nearly scale invariant large-scale power spectrum 
from inflation. The cross term (5i4>82(j)} is of order ((Si<p) 4 } and represents a negligible "loop" correction to the scale- 
invariant spectrum from inflation. (The cross term does not involve the iso-inflaton since 5\4> and S\x arc statistically 
independent.) Finally, there is a contribution {(52<t>) 2 ) which involves terms of order (x 4 ) coming from rescattering 
and terms of order ((5i0) 4 ) which represent (more) loop corrections to the scale- invariant spectrum from inflation. 
Thus, we can schematically write 

pm = pr{k) [i + (bops)] + p^ c (k) 

Here P™ ~ fc" s_1 is the usual nearly scale invariant spectrum from inflation and P^ csc is the bump-like contribution 
from rescattering and IR cascading which we have studied in t he pr evious section. The "loop" corrections to P™ c (fc) 
have been studied in detail in the literature (see, for example, [103]) and are known to be negligible in most models. 

We can also make a similar schematic decomposition of the bispectrum by considering the structure of the 3-point 
correlator ((8<fi) 3 ). Following our previous line of reasoning, it is clear that the dominant contribution comes from 
rescattering and is of order (x 6 )- The terms involving (((5i</>) 3 ), on the other hand, represent the usual nongaussianity 
generated during single field slow roll inflation and are known to be small [§8[ . 

D. The Curvature Perturbation 

Ultimately one wishes to compute not the field perturbations 8 n <p, S n x, but rather the gauge invariant curvature 
fluctuation, £. We expand this in perturbation theory in the usual manner 

C = Ci + \(2 (99) 

In jlOOl ] Malik has derived expressions for the large scale curvature perturbation in terms of the Sasaki-Mukhanov 
variables at both first and second order in perturbation theory. We remind the reader that in the flat slicing (which 
we employ) the Sasaki-Mukhanov variable for each field simply coincides with the field perturbation (i.e. - = 6(f) 
and Q x = S\)- 

At first order in perturbation theory the iso-inflaton does not contribute to the curvature perturbation (since 
(x) = 0) and we have 



(100) 
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At seco nd o rder in perturbation theory the expression for the curvature perturbation is more involved. Using the 
results of [l00j | and working to leading order in slow roll parameters we find 



C2 



H 

> 
1 



3(</>') 2 
1 



m 

.^[(SiW + aWfrtf] (101) 

Let us discuss the various contributions to this equation. The third line contributes to the nongaussianity of the 
vacuum fluctuations during inflation. These terms are known to be negligible [l0l - fl2l l98j and, indeed, one may 
explicitly verify that (I101[) would predict Jnl ~ 0(e,rj) in the absence of particle production. 

Next, we consider the second line of (|101[) . This represents the direct contribution of the gaussian fluctuations 
8\x to the curvature perturbation. This contribution is tiny since the x particles are extremely massive for nearly 
the entire duration of inflation and hence Six ~ a~ 3 / 2 (see also (3(| for a related discussion). The smallness of this 
contribution to £ can be understood physically by noting that the super-horizon iso-curvature fluctuations in our 
model are negligible. 

Finally, let us consider the contribution on the first line of (I101[) . This contribution is the most interesting. To 
make contact with observations we must compute the curvature perturbation at late times and on large scales. In 
section HVl we have already shown that 5 n <j> is constant on large scales and at late times for both n = 1 and n = 2. 
This is the exp ecte d result: the curvature fluctuations are frozen far outside the horizon and in the absence of entropy 
perturbations 106]. 10 Hence S2<t>' is completely negligible and the first term on the first line of (|101l) must dominate 
over the second term. We conclude that ,at late times and on large scales, the second order curvature perturbation is 
very well approximated by 

C 2 = -^20+--' (102) 

In summary, we have shown that the power spectrum of curvature fluctuations from inflation in the model ([5]) is 
trivially related to the power spectrum of inflaton fluctuations 

P c (fc) <* = 2^P*{k) (103) 

at both first and second order in cosmological perturbation theory. This relation is valid at late times and for scales 
far outside the horizon. The curvature spectrum (|103[) may be written as 

P c (k) = P c vac (fc) [1 + (loops)] + Pf sc (fc) (104) 

The power spectrum of the inflaton vacuum fluctuations agrees with the usual result obtained in linear theory (94j 

H 2 ( k ^ 2r, ~ 6e 



In (|104[) we have schematically labeled the corrections arising from the second line of (I101[) and the source as 
"loop" . These are nongaussian corrections to the inflaton vacuum fluctuations arising from self-interactions of the 
inflaton and also the nonlinearity of gravity. Such corrections are negligible. The most interesting contribution to the 
power spectrum (|104j) is due to rescattering, P£ esc (k). This quantity is proportional to our previous result (|71j) . 

In passing, notice that the bispectrum (defined by (174^ *) of inflaton fluctuations will differ from the bispectrum 
B of the curvature fluctuations (defined by (|3])) only by a simple re-scaling: 

H\ 3 „ „ , 1 



Bi l; . - ( -r ) B+(ki) = ~ {2e)mM f *(h) (106) 



The dominant contribution to B^ comes from rescattering effects (S2<j> 3 ) ~ (Six 6 )- 
The analysis of this section justifies our neglect of metric fluctuations in section HV 



Note that, in some cases, the curvature fluctuations may evolve significantly after horizon exit [l07l . Il0§ |. This is a concern in models 
where there are significant violations of slow-roll. In section [HI CI we have already shown that the transient violation of slow roll has a 
negligible effect on the curvature fluctuations in our model. Hence, the result £ n ~ 5 n <j> ~ const far outside the horizon is consistent 
with previous studies. 
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VI. OBSERVATIONAL CONSTRAINTS ON PARTICLE PRODUCTION DURING INFLATION 

In the previous sections of this review, we have demonstrated that particle production and IR cascading in the 
model Q leads to a bump-like feature in the primordial power spectrum. We now aim to determine the extent 
to which such a spectral distortion is compatible with current cosmological data. The results of this section first 
appearing in the paper Q. 

One motivation for this study is to determine whether features generated by particle production during inflation 
can explain some of the anomalies in the observed primordial power spectrum, P(k). A number of different studies 
have hinted at the possible presence of some localized features in the power spectrum [H, [zIHIIj], which arc 
not compatible with the simplest power law P(k) ~ fc™ 3 ^ 1 model. Although such glitches may simply be statistical 
anomalies [871 ) , there is also the tantalizing possibility that they represent a signature of primordial physics beyond the 
simplest slow roll inflation scenario. Forthcoming polarization data may play a crucial role in distinguishing between 
these possibilities (54j . However, in the meantime, it is interesting to determine the extent to which such features 
may be explained by microscopically realistic inflation models, such as ([5])- 

To answer this question we provide a simple semi-analytic fitting function that accurately captures the shape of 
the feature generated by particle production and IR cascading during inflation. Next, we confront this modified 
power spectrum with a variety of observational data sets. We find no evidence for a detection, however, we note that 
observations are consistent with relatively large spectral distortions of the type predicted by the model (|4|). If the 
feature is located on scales relevant for Cosmic Microwave Background (CMB) experiments then its amplitude may 
be as large as 0(10%) of the usual scale-invariant fluctuations, corresponding to g 2 ~ 0.01. Our results translate 
into a 0o-dependent bound on the coupling g 2 which is crucial in order to determine whether the nongaussian signal 
associated with particle production and IR cascading is detectable in future missions. 

We also consider the more complicated features which result from multiple bursts of particle production and IR 
cascading. Such features are a prediction of a number of string theory inflation models, including brane/axion 
monodromy (58l - [60l |. For appropriate choice of the spacing between the features, we find that the constraint on g 2 in 
this scenario is even weaker than the single-bump case. 

A. A Simple Parametrization of the Power Spectrum 

In Q it was shown that particle production and IR cascading during inflation in the model Q generates a bump-like 
contribution to the primordial power spectrum. As shown in Fig. [51 this feature can be fit with a very simple function 
-Ptmmp ~ fc 3 e~ 7rfc /( 2fc *). The bump- like contribution from IR cascading is complimentary to the usual (nearly) scale- 
invariant contribution to the primordial power spectrum P vac ~ fc™ s_1 coming from the quantum vacuum fluctuations 
of the inflaton. The total, observable, power spectrum in the model Q is simply the superposition of these two 
contributions: P(k) ~ k ns ~ x + k 3 e~ vk ^ 2k *\ (See equation (|65p.) This simple formula can be motivated from 
analytical considerations @ and provides a good fit to lattice field theory simulations near the peak of the feature 
and also in the IR tail. 11 

It is straightforward to generalize this discussion to allow for multiple bursts of particle production during inflation. 
Suppose there are multiple points (p = (pi (i = 1, • • • , n) along the inflationary trajectory where new degrees of freedom 
Xi become massless: 

" 2 

4nt = -^f (107) 

i=0 

For each instant ti when <fi — cf>i there will be an associated burst of particle production and subsequent rescattering 
of the produced massive \i on? the condensate <fi(t). Each of these events proceeds as described above and leads to a 
new bump-like contribution to the power spectrum. These features simply superpose owing to that fact that each field 
Xi is statistically independent (so that the cross terms involving XiXj with i ^ j in the computation of the two-point 
function must vanish). Thus, we arrive at the following parametrization of the primordial power spectrum in models 



11 This fitting formula docs not capture the small oscillatory structure in the UV tail of the feature (see Fig.QJJl which does not concern us 
since that region is not phenomenologically interesting. 



30 



with particle production during inflation: 

where A s is the amplitude of the usual nearly scale inva riant vacuum fluctuations from inflation and ko is the pivot, 
which we choose to be ko = 0.002 Mpc -1 following [l09j |. The constants A, depend on the couplings g 2 and measure 
the size of the features from particle production. We have normalized these amplitudes so that the power in the i-th 
bump, measured at the peak of the feature, is given by Ai. The location of each feature, hi, is related to the number 
of e-foldings N from the end of inflation to the time when the i-th burst of particle production occurs: roughly 
hi(ki/H) ~ Ni where N = Ni at the moment when (j> = fa. From a purely phenomenological perspective the locations 
ki are completely arbitrary. 




FIG. 9: The bump-like features generated by IR cascading. We plot the feature power spectrum Pbump(fc) obtained from fully 
nonlinear lattice field theory simulations (the red points) and also the result of an analytical calculation (the dashed blue curve) 
using the formalism described in §]. We also superpose the fitting function ~ k 3 e~ nk /( 2fc *' (the solid black curve) to illustrate 
the accuracy of this simple formula. 



We compare (|I08|) to lattice field theory simulations in order to determine the amplitude Ai in terms of model 
parameters. We find 

A 4 -10- 6 5 f /4 (109) 

The power-law form of this relation was determined by inspection of numerical results, however, it can also be 
motivated by analytical considerations. Assuming standard chaotic inflation V = m 2 (j) 2 /2 we have tested this formula 
for g 2 — 1,0.1,0.01, taking both 0o = 2V&ttM p and 4>o = 3.2^/8ttM p . We found agreement up to factors order unity 
in all cases. 

Theoretical consistency of our calculation of the shape of the feature bounds the coupling as 10~ 7 ^ g 2 ^, 1 jg|. 
Hence, the power spectrum (|108[) can be obtained from sensible microphysics only when 10~ 20 <, Ai <, 10~ 6 . This 
constraint still allows for a huge range of observational possibilities: near the upper bound the feature is considerably 
larger than the vacuum fluctuations while near the lower bound the feature is completely undetectable. 

Note that for each bump in (| 108[) the IR tail Pbump — > k 3 as k — > is similar to the feature considered by Hoi, 
Cline & Holder in [84| , consistent with causality arguments about the generation of curvature perturbations by local 
physics. 



B. Data Sets and Analysis 



The primordial power spectrum for our model is parametrized as (|108l) . Our aim is to derive observational con- 
straints on the various model parameters A s , n s , ki and Ai using CMB, galaxy power spectrum and weak lensing 
data. To this end we use the cosmoMC package |l 1 0i | to run Markov Chain Monte Carlo (MCMC) calculations to 
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determine the likelihood of the cosmological parameters, including our new parameters Ai and fcj. We employ the 
following data sets. 

Cosmic Microwave Background (CMB) 

Our complete CMB data sets include WMAP-5yr fiol [ml. B OOMERANG [mWlllj . ACBAR [TTsMTTsjj . 
CBI [nHHl , VSA [123, DASI [Hi IHl. and MAXIMA~]l23. We have included the Sunyaev-Zeldovic (SZ) 



secondary anisotropy |127lll28ll f or WMAP-5yr, ACBAR and CBI data sets. The SZ template is obtained from 
hydrodynamical simulation [129]. Also included for theoretical calculation of CMB power spectra is the CMB 
lensing contribution. 

Type la Supernova (SN) 

We employ the Union Supernova la data (307 SN la samples) from The Supernova Cosmology Project |l30j . 

Large Scale Structure (LSS) 

The 2dF Galaxy Redshift Survey (2dFGRS) data [l3l| and Sloan Digital Sky Survey (SDSS) Luminous Red 



Galaxy (LRG) data release 4 [132j are utilized. 

Note that we have used the likelihood code based on the non-linear modeling by Tegmark et al. |l32j (marginal- 
izing the bias b and the Q parameter). However with a large bump in the linear power spectrum, this naive 
treatment may be not sufficient to characterize the non-linear response to the feature on small scales. Ideally, 
this should be obtained from N-body simulations, however, such a study is beyond the scope of the current 
work. 

There are several other caveats on our results in the high-fc regime. First, we assume linear bias for the galaxies, 
which may not be entirely safe at sufficiently small scales. Moreover, sharp features in the matter power spectrum 
can cause sharp features in the bias as a function of k. 

Keeping in mind these caveats our constraints on small scales k >, 0.1 Mpc -1 should be taken with a grain of 
salt and considered as accurate only up to factors order unity. 

Weak Lensing (WL) 

Five WL data sets are used in this paper. The effective survey area A e f{ and galaxy number density n e ff of each 
survey are listed in Table HI 



TABLE I: Weak Lensing Data Sets 



Data sets 


Aeff 


n e g 




(deg 2 ) 


(arcmin -2 ) 


COSMOS [133, 1391 


1.6 


40 


CFHTLS-wide f!34. 1351 


22 


12 


GaBODS [136, .137] 


13 


12.5 


RCS [136. 1371 ~ 


53 


8 


VIRMOS-DESCART [135. 1381 


8.5 


15 



For COSMOS data we use the CosmoMC plug-in written by Julien Lesgourgues [l39j . modified to do numerical 
marginalization on three nuisance parameters in the original code. 

For the other four weak lensing data sets we use the likelihood given by [140} . To calculate the likelihood we 
have written a CosmoMC plug-in code, with simplified marginalizatio n on the parameters of galaxy number 
density function n(z). More details about this plug- in can be found in [l4l| . 

As for the LSS data, for small scales k ^ 0.1 Mpc" 1 there is the caveat that the nonlinear evolution of the power 
spectrum in the presence of bump-like distortions may not be treated accurately. 



C. Observational Constraints: A Single Burst of Particle Production 



We now present our results for the observational constraints on particle production during inflation, assuming two 
different scenarios. 
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Nonlinear scales 
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FIG. 10: Marginalized posterior likelihood contours for the parameters Air and km (the magnitude and position of the feature, 
respectively) in the single-bump model. Black and grey regions correspond to parameter values allowed at 95.4% and 99.7% 
confidence levels, respectively. At small scales, to the right of the dashed vertical line, our results should be taken with a 
grain of salt since the nonlinear evolution of the power spectrum may not be modeled correctly in the presence of bump-like 
distortions. 



TABLE II: Constraints on the standard ("vanilla") cosmological parameters for the single-bump model. All errors are 95.4% 
confidence level 
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The minimal scenario to consider is a single burst of particle production during inflation, which corresponds to 
taking n = 1 in ()10T[) . The power spectrum is given by (|108[) with n = 1 and, with some abuse of notation, 
we denote fci = &ir and A\ = Air. The prior we have used for Air is < Air < 25 x 10~ 10 , and for fcrR is 
—9.5 < ln[A:/Mpc _1 ] < 1. The former condition ensures that the bump-like feature from IR cascading does not 
dominate over the observed scale invariant fluctuations while the latter is necessary in order to have the feature in the 
observable range of scales. In Fig. [10] we plot the marginalized posterior likelihood for the new parameters Air and 
A:ir describing the magnitude and location of the bump while in Table HI] we give the best fit values for the remaining 
(vanilla) cosmological parameters. 

For very large scales ;$ Gpc -1 , the data do not contain much information (due to cosmic variance) and hence the 
constraint on any modification of the power spectrum is weak. In this region the spectral distortion may be larger than 
100% of the usual scale invariant fluctuations and couplings g 2 order unity are allowed. For smaller scales k Gpc -1 
the constraints are stronger and we have, very roughly, Air/A s < OA corresponding to g 2 < 0.01. For very small 
scales, k > 0.1 Mpc" 1 our constraints should be taken with a grain of salt since the nonlinear evolution of the power 
spectrum may not be modeled correctly in the presenc e of bump-like distortions. At small scales nonlinear effects 
tend to wipe out features of this type (see, for example, [l42j |) and hence observational constraints for k ^ 0.1 Mpc -1 
may be weaker than what is presented in Fig. 1101 Note that in most of this nonlinear regime we find essentially no 
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FIG. 11: Marginalized posterior likelihood contours for the parameters Air and A (the feature amplitude and spacing, respec- 
tively) of the multiple-bump model. Black and grey regions correspond to values allowed at 95.4% and 99.7% confidence levels, 
respectively. 



constraint on Air, which is consistent with what would be expected in a more comprehensive treatment. 

The IR cascading bump in the primordial power spectrum will be accompanied by a corresponding nongaussian 
feature in the bispectrum, that will be discussed in more detail in the next section. From the perspective of potentially 
observing this signal it is most interesting if this feature is located on scales probed by CMB experiments. (There is 
also th e fascinat ing possibility that the nongaussianity from IR cascading could show up in the large scale structure as 
in (2ll Il43l4l45j |. We leave a detailed discussion to future studies.) To get some intuition into what kinds of features 
in the CMB scales are still allowed by the data we focus on an example with Air — 2.5 x 1CP 10 which, using (I109p . 
corresponds to a reasonable coupling value g 2 ~ 0.01. We take the bump to be located at fciR = 0.01 Mpc -1 and fix 
the remaining model parameters to A s = 2.44 x 10~ 9 , n s = 0.97 (which are compatible with the usual values). This 
sample bump in the power spectrum is illustrated in the top panel of Fig. [3] and is consistent with the data at 2a. In 
the bottom panel of Fig. [3] we plot the associated angular CMB TT spectrum. This example represents a surprisingly 
large spectral distortion: the total power in the feature as compared to the scale invariant vacuum fluctuations is 
-Pt>ump/-fvac ~ 0.1, evaluated at the peak of the bump. 

D. Observational Constraints: Multiple Bursts of Particle Production 

Next, we consider a slightly more complicated scenario: multiple bursts of particle production leading many localized 
features in the power spectrum. For simplicity we assume that all bumps have the same magnitude Aj = Am and 
we further suppose a fixed number of e-foldings SN between each burst of particle production. This implies that 
the location of the i-th bump will be given by k% — e^~^'^ki where A is a model parameter controlling the density 
of features. We take the number of bursts, n, to be sufficiently large that the series of features extends over the 
whole observable range. In the next section we will see that these assumptions are not restrictive and that many 
well-motivated models are consistent with this simple set-up. 

Our multi-bump model, then, has three parameters: Am, k\ and A. We take the prior on the amplitude to be 
Air < 25 x 10~ 10 as in section fVI CI If the features are very widely spaced, A >, 1, then the constraint on each 
bump will obviously be identical to the results for the single-bump case presented in the section IVI CI Hence the 
most interesting case to consider is A < 1 so that the bumps are partially overlapping. Our prior for the density of 
features is therefore < A < 1. Finally, the location of the first bump will be a historical accident in realistic models, 
hence we marginalize over all possible values of k\ and present our constraints and 2-d likelihood plots in the space 
of Air and A. This marginalized likelihood plot is presented in Fig. QTJ In table Hm we present the best-fit values for 
the vanilla cosmological parameters. 

From the likelihood plot, Fig. [TTJ there is evidently a preferred value of the feature spacing, roughly A ~ 0.75, for 
which the constraints are weakest. This can be understood as follows. For very high density A — > the localized 
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TABLE III: constraints on the standard ("vanilla") cosmological parameters for the multiple-bump model. All error bars are 
95.4% confidence level. 
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features from IR cascading smear together and the total power spectrum (|108p is P(k) ~ A s (k/ko) ns ^ 1 + C where the 
size of the constant deformation scales linearly with the density of features: C oc A" 1 . Therefore, the upper bound on 
the amplitude should scale linearly with A. Indeed, this linear trend is very evident from Fig. [11] in the small- A 
regime. This linear behaviour must break down at some point since as the features become infinitely widely spaced 
the constraint on Am must go to zero. This explains the bump in the likelihood plot, Fig. [TTJ near A ~ 0.75. 

In passing, notice that the behaviour P(k) ~ A s {k/ko) ns ~ 1 + C for A <§C 1 also explains why the best-fit A s in 
table IIIII is somewhat lower than the standard value and why the spectral tilt n s — 1 is somewhat more red. 

To get some intuition for the kinds of multi-bump distortions that are allowed by the data, we consider an example 
with A\n = 1 x 10~ 9 , A = 0.75 and fix the vanilla parameters to A s = 1.04 x 10 -9 , n s = 0.93. This choice of 
parameters is consistent with the data at 2a and corresponds to a reasonable coupling g 2 ~ 0.02. In Fig. Q2] we plot 
the primordial power spectrum P(k) and also the CMB TT angular power spectrum for this example. 

VII. PARTICLE PHYSICS MODELS 

In previous sections of this review we have studied in some detail the observational signatures of inflationary particle 
production in the model ^ and we have also derived observational constraints on the model parameters g 2 and <f>Q. 
We now aim to provide some explicit microscopic derivations of the model ([5]) from popular theories of particle physics 
including string theory and SUSY. Our examples are meant to be illustrative, rather than exhaustive. From the low 
energy perspective interactions of the type (j4|) are completely generic and hence we expect that many microscopic 
derivations may be possible. 

A. Open String Inflation Models 

String theory inflation models may be divided into two classes depending on the origin of the inflaton: closed string 
models and open string models. In the former case the inflaton is typically a geometrical modulus assoc iated with 
the compactification manifold (e xam ples include racetrack inflation [l46j ] . Kahler modulus inflat ion Jl47l| . Roulette 
inflation [l48j | and fibre inflation |l49l |). Reheating after closed string inflation has been studied in 150]. On the other 
hand, in the case of open string models the inflaton is typically the posit ion modulus of some mobile D-brane 12 livin g 
in the compactifi catio n manifold (examples include brane inflation [l5ll | such as the warped KKLMMT mod el [152| | , 
D3/D7 inflation jl53| and DBI inflation [34| ). Reheating after open string inflation has been studied in |l54j . 

In open string inflation models there may be, in addition to the mobile inflationary brane, some additional "specta- 
tor" branes. If the mobile brane collides with any spectator brane during inflation then some of the stretched string 
states between these branes will become massless at the moment when the branes are coincident [36|, [6l[ , precisely 
mimicking the interaction @. Thus, we expect particle production, IR cascading and the bump- like features described 
above to be a reasonably generic prediction of open string inflation. 



12 One notable exception is inflation driven by the open string tachyon, for example nonlocal string field theory models [29.1321 15^1 ■ 
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FIG. 12: The top panel shows a sample multiple-bump distortion with amplitude Am = 1 x 10 which corresponds to 
g 2 ~ 0.02. The feature spacing is A = 0.75. We take the vanilla parameters to be A s = 1.04 x 10 -9 , n a — 0.93 so that the scale 
of inflation is slightly lower than in the standard scenario and the spectral tilt is slightly redder. The bottom panel shows the 
CMB angular TT power spectrum for this example. 



B. String Monodromy Models 



A concrete example of the heuristic scenario discussed in the last subsection is provided by the brane monodromy 
and axion monodromy string theory inflation models proposed in [58-63]. In the original brane monodromy model 
[5c| one considers type IIA string theory compactified on a nil manifold that is the product of two twisted tori. The 
metric on each of these twisted tori has the form 

ds 2 

— = L 2 Ui dul + L 2 U2 duj + L 2 x (dx' + Mu x du 2 f (110) 

where x' — x — M-U1U2 and M is an integer flux number. The dimensionless constants L Ul , L U2 and L x determine 
the size of the compactification. 

Inflation is realized by the motion of a D4-brane along the direction U\ of the internal manifold. The D4 spans our 
large 3-dimensions and wraps a 1-cycle along the direction U2 of the internal space. The size of this 1-cycle, in string 
units, is given by 




Hence, the brane prefers to minimize its world- volume by moving to the location u\ = where this 1-cycle has minimal 
size. This preference gives a potential to the D4-brane position which goes like V ~ u\ in the large u± regime that is 
relevant for large field inflation. 

In [36j it was shown that this scenario allows for the inclusion of a number of spectator branes stabilized at positions 
u\ = i/M (with i integer) along the inflationary trajectory. As the mobile inflationary D4 rolls through these points 
particles (strings) will be produced and the resulting distribution of features will look nearly identical to the simple 
multi-bump scenario studied in section IVI Dl To see this, let us now determine the distribution of features that is 
predicted from brane monodromy inflation. The canonical inflaton <p can be related to the position of the mobile D4 
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as 

ct> = Bu\ /p (112) 

where B is a constant with dimensions of mass that depends on model parameters. Hence, the effective potential 
during inflation has the power-law form 

V(cP) = M 4 -P0P (113) 

For the simplest scenario described above one has p = 2/3. However, the formulas (| 1 1211 13)) still hold for the variant 
considered in [58j with p = 2/5 as long as one replaces u± by a more complicated linear combination of coordinates. 
These relations also hold for axion monodromy models [591 ] with p — 1 and u\ replaced by the axion, c, arising from a 
2-form RR potential C 1 - 2 ' integrated over a 2-cycle E 2 . For all models of the form (I113[) the number of e- foldings N 
from (f> = (j)(N) to the end of inflation (which occurs at <f> ~ p/V% when the slow roll parameter e((f>) = 1) is given by 

= 1 2 (AQ p 
2p Ml 4 

2pM 2 1 4 V ' 

Since the spectator branes are located at u\ — i/M the bursts of particle production must occur at times N = iVj 
where 

The location k — ki of the i-th feature is defined, roughly, by the scale leaving the horizon at the moment N = A/j. 
Hence, the distribution of features for brane/axion monodromy models is given by 



In 



B 2 i 2 'v - V - (116) 
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with p = 2/3 or p = 2/5 for brane monodromy and p = 1 for axion monodromy. In (|116[) the dimensionless number 
B depends on model parameters. 

Although the distribution of features (|116p is not exactly the same as the evenly space distribution considered 
subsection IVI D[ the two are essentially indistinguishable over the range of scales which are probed observationally 
(corresponding to roughly 10 e-foldings of inflation). The reason for this is simple: the inflaton is nearly constant 
during the first 10 e-foldings of inflation and hence 5N ~ 5<j) ~ 5u\ within the observable region. It follows that 
ki = e( i_1 ) A fci to very good approximation for a huge class of models. To see this more concretely in the case at 
hand, let us compute dN/dui from (|114[) . It is straightforward to show that 



dN 1 1 ( B 



dux P p [2t{ ( j ) )] l ' p/2 \M P 

where 



(117) 



2 V V J 2 \ 

is the usual slow roll parameter. Observational constraints on the running of the spectral index imply that e(</>) cannot 
change much over the observable 10 e-foldings of inflation. Since dN/du\ = const to very high accuracy it follows 
trivially that N — N(u\) is very close to linear and ki = e^ _1 ^ A fci as desired. 

In the context of axion monodromy inflation models [59j the multiple bump features discussed here will be com- 
plimentary to the oscillatory features described in [60l ] which result from the sinusoidal modulation of the inflaton 
potential by instanton effects. If the bursts of particle production are sufficiently densely spaced, then signal from IR 
cascading may appear oscillatory, however, it differs from the effect discussed in [60] in both physical origin and also 
in functional form. 

Let us now estimate the effective value of the couplings gf appearing in the prototype interaction (|107l) that are 
predicted from the simplest brane monodromy model. A complete calculation would involve dimensionally reducing 
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the DBI action describing the brane motion and requires knowledge of the full 10-dimensional geometry with the 
various embedded branes. For our purposes, however, a simple heuristic estimate for the collision of two D4-branes 
will suffice. When N D-branes become coincident the symmetry is enhanced from U(1) N to a U(N) Yang Mill gauge 
theory The gauge coupling for this Yang Mills theory is given by 

9ym = y ~\^ (H9) 

where L is the volume of the 1-cycle that the D4 branes wrap and is given by ([5]). If the inflationary brane is at 
position Mi and the i-th spectator brane is at u\ t i then the distance between the two branes is given by 

d 2 = a' L 2 Ui ( Ul - Ul ,i) 2 (120) 

The mass of the gauge bosons corresponding to the enhanced symmetry is 

Ml = .vmt^TTW = *f>*-""> 2 (121) 



(27r) 2 (a') 2 



>JL* +L%M*ui 



To put this in the prototype form M 2 — g 2 ((f> — 4>i) 2 we must first convert to the canonical variable <j> using the formula 
(TTT2|) with p = 2/3 and 

M 1 / 2 L j}J 2 

B = j. 2 J= r 122 

Next, we must Taylor expand the resulting equation about the minimum <f> — <fn. We find 

M? = g 2 {4> - cf> t ) 2 + ■ ■ ■ (123) 
16<7 2 tt 4 1 



ML x U\.i ^ L 2 a+L 2 M 2 U 2 

16g 2 7r 4 1 
L x i y/Ll 2 +L*i* 



(124) 



where on the second line of (|124p we have used the fact that u\.i — i/M (with i integer) in the simplest models. We 
see that the effective couplings g 2 become larger as the D4 unwinds during inflation. (The apparent divergence for 
u\ : i = in the formula (I124p is an artifact of the fact that the relation (|112[) is not valid at small values of U\, This 
will not concern us here since inflation has already terminated at the point that our formulas break down.) 

To compute the amplitude of the bump-like feature produced by brane monodromy inflation we should take into 
account also combinatorial factors. When two branes become coincident the symmetry is enhanced from U(l) 2 to 
U{2) so there are 2 2 — 2 = 2 additional masslcss spin-1 fields appearing at the brane collision. Thus, using equation 
(|109p , the amplitude of the feature that will be imprinted in the CMB is 

A i>eS = 2 x (2 2 - 2) x [lCT 6 • g\ 5/4 ] (125) 

where the extra factor of 2 counts the polarizations of the massless spin-1 fields. This combinatorial enhancement can 
be much larger if the inflationary brane collides with a stack of spectators. 

The above discussion is predicated on the assumption that the original brane monodromy set-up 1 is supplemented 
by additional spectator branes. This may seem like an unnecessary contrivance, however, in order for this model 
to reheat su ccess fully it may be necessary to include spectator branes. For example, with the reheating mechanism 
proposed in [155| semi-realistic particle phenomenology can be obtained by confining the standard model (SM) to a D6 
brane which wraps the compact space. In order to reheat into SM degrees of freedom we orient this brane so that its 
world- volume is parallel to the mobile (inflationary) D4. In this case the end of inflation involves multiple oscillations 
of the D4 about the minimum of its potential. At each oscillation the D4 collides with the D6 and SM particles are 
produced by parametric resonance preheating (68| . However, due to the periodic structure of the compactification, 
D4/D6 collisions will necessarily occur also during inflatiomjeading to IR cascading features in the CMB. 

The timing of these D4/D6 collisions was com put ed in [1551 ] for the minimal p — 2/3 brane monodromy model, 
assuming the same choices of parameters used in [581 ]. For this particular case there is only one collision (and hence 
one feature) during the first 10 e-foldings of inflation and the phenomenology is essentially the same as that considered 
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in subsection IVI CI What is the amplitude of this f eatur e? Assuming, again, the parameters employed in [58| and 
noting that the first collision takes place at i = 13 [155| equation f|l 24[) gives g\ = 0.001. From (|125[) we find the 
effective amplitude of the feature to be Ax ie g/A s = 0.01. This value is well within the observational bounds derived 
in subsection I VI CI however, it should be within reach of future missions (88j . 

We stress that the conclusions in the previous paragraph apply only for the particular choice of model parameters 
employed in |58j . There exist other consistent parameter choices for which the simplest brane monodromy model 
predicts a much higher density of features with much larger amplitude. 

Note that both brane and axion monodromy models may be used to realize trapped inflation [36J. Here we are 
restricting ourselves to the large-field regime where the potential V = fi 4 ~ p cf)' p is flat enough to drive inflation without 
the need for trapping effects. For a given choice of parameters one should verify that this classical potential dominates 
over the quantum corrections from particle production. 



C. A Supersymmetric Model 



Another microscopic realization of multiple bursts of particle production and IR cascading during inflation which 
does not rely on string theory can be obtained from the so-called "distributed mass" model derived in [57j with warm 
inflation [63j in mind, however, the theory works equally well for our scenario. This model is based on Af = 1 global 
SUSY and allows for the inclusion of multiple points along the inflationary trajectory where both scalar degrees of 
freedom and also their associated fermion super-partners become massless. The distribution of features in this set-up 
is essentially arbitrary. 



VIII. NONGAUSSIANITY FROM PARTICLE PRODUCTION AND IR CASCADING 



We have seen that the dynamics of particle production, rescattering and IR cascading during inflation in the model 
([5]) leads to a bump-like contribution to the primordial power spectrum of the observable cosmological fluctuations. 
However, this dynamics must also have a nontrivial impact on nongaussian statistics, such as the bispectrum. Indeed, 
it is evident from the analytical formalism presented in sections IIVI and fVl that the inflaton fluctuations S<fi generated 
via the rescattering diagram Fig. [T] may be significantly nongaussian. (This is evident since the particular solution of 
equation (f5Tj) is bi-linear in the Gaussian field y.) In this section we characterize this nongaussianity. The results of 
this section appeared for the first time in [1, d, 0| . 



A. Probability Density Function and Cummulants 



As discussed in the introduction, nongaussianity is often characterized by the size, shape and running of the 
bispectrum B(k{). A popular measure of the size of the nongaussianity is the magnitude of B(ki) on equilateral 
triangles. This approach has the advantage of being straightforward, however, it may give misleading results when 
one wishes to compare bispectra with different shapes 20] or scaling behaviour. This is particularly true in models, 
such as ([5]) , where the bispectrum is large only for triangles with a characteristic size and where the nongaussian part 
of £ is uncorrelated with the gaussian part. A more meaningful single number to compare between models is the 
skewness (defined below) , which contains information ab out t he bispectrum integrated over all scales (up to some UV 
smoothing scale) and all shape configurations. (See also [l58j | for a related discussion and alternative methodology.) 

In what follows it will be useful to define a normalized curvature perturbation 

£=— (126) 
where £ is the usual gauge invariant curvature perturbation and is the variance 

A = '■(') = J (££,(£^ <o-.a.> < 127 » 

If C is generated by the quantum vacuum fluctuations of the inflaton field then we have o~q = 10~ 9 / 2 . 

Let us define the probability density function (PDF), P(£), as the probability that the (normalized) curvature 
perturbation has fluctuations of size £. The n-th central moment of the PDF is simply 

/+oo 
CP(Z)d£ (128) 
-oo 
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The n-th dimensionless cummulant k n is the connected n-point function. For (£) = the first few non-vanishing 
cummulants are: 

« 3 = (C 3 > (129) 

« 4 = (e 4 > - 4(e 2 > 2 (iso) 

For a gaussian PDF we have k n = for n = 3, 4, • • • , hence the cummulants provide a measure of nongaussianity. 
In the case that ( (and hence £) is a free field then the central limit theorem implies that the PDF is gaussian 

P(£K = -g=e-^ (131) 

V Z7T 

This expression admits a simple generalization to the case where P(£) is close to (but not exactly) gaussian. This 
generalization is the Edgeworth expansion 



V Z7T 



l + yP 3 (C) 



(132) 



(See [21| for a review and [1591 ] for an alternative derivation.) In (|132l) H n (Q denotes the Hermite polynomials of order 
n. Note that the result (|132[) is an expansion in cummulants. This expression provides an accurate approximation of 
the true PDF provided the cummulants are well ordered in the sense that 

1 > \k 3 \ > > ••• (133) 

See [HI 164] for more discussion on the structure of the correlation functions and the validity of cosmological 
perturbation theory. 

Finally, it remains to relate the cummulants to the correlation functions (such as the bispectrum) of the primordial 
cur vatur e perturbation. We can write the dimensionless skewness as an integral over the 3-point correlation function 

fHfTH 



1 f d 3 h d 3 k 2 d 3 k 3 

~ ^ J (2^)3/2 ( 27r )3/2 (2^)3/2 (Ck 1 Ck 2 Ck 3 )c (134) 



A similar expression may be derived for the dimensionless kurtosis, S4. From (|134|) it is clear that the skewness 
provides a measure of the magnitude of the bispectrum integrated over all triangles. Hence, this provides a reasonable 
single number to compare the size of nongaussianity between different models. Moreover, the skewne ss (and highe r 
cum mulants) are of observational interest since these determine the probability of rare fluctuations [2l|, ImMIM 1169- 
162]. Observables such as the abundance of collapsed objects may therefore be used to constrain k 3 , K4, etc. 

We have reviewed the derivation of the Edgeworth expansion (|132|) only to illustrate that the dimensionless skewness 
(11341) encodes information about the bispectrum integrated over all momenta. However, our actual calculation of the 
PDF in the next subsection will be fully non-perturbative and hence does not rely on the validity of (|132p or (|133|) . 



B. Magnitude of the Nongaussianity from IR Cascading 

We have argued in the last subsection that the moments of the PDF provide a useful tool for quantifying the size of 
nongaussianity and drawing comparisons between microscopic models whose bispectra may have very different shapes 
and scaling properties. Our goal now is to compute the skewness and kurtosis generated by particle production, 
rescattering and IR cascading in the model ([5]). We will then compare this to the skewness that would be generated 
by more familiar constructions, such as the local model ( = C s + f/wLCg- 

We extract the PDF of d(f> from our HLattice simulations by measuring the fraction of the simulation box which 
contains the fluctuation field 5<f> at a particular value. Notice that this approach is completely nonperturbative: it 
does not rely on the validity of the Edgeworth expansion, nor does it assume anything about the size or ordering of 
the cummulants. This procedure implicitly puts a IR cut-off at the box size L and a UV cut-off at the lattice spacing, 
A -1 . Since the nongaussian effects in our model are strongly localized in Fourier space, our quantitative results are 
largely insensitive to L and A. 

In Fig. [13] we plot our numerical result for the PDF of the inflaton fluctuations generated by rescattering and 
IR cascading. In order to make the physics of inflationary particle production clear, we have subtracted off the 
contribution coming from the usual vacuum fluctuations of the inflaton. That is, the PDF in Fig. Q2] is associated 
only with the contribution 6<j) iesc in equation (1171) . 
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FIG. 13: The PDF of the inflaton fluctuations generated by rescattering and IR cascading, at a series of different values of the 
scale factor, a. The dotted black curve shows a Gaussian fit at late times and we have normalized the scale factor so that a — 1 
at the moment when particle production occurs. For illustration, we have chosen g 2 = 0.1 and a standard chaotic inflation 
potential V(cj>) = rn 2 cj> 2 /2. 



We can understand physically the behaviour of PDF plotted in Fig. [13] Shortly after the initial burst of particle 
production the inflaton perturbations 5(j) are extremely nongaussian, due to the sudden appearance of the source term 
J oc x 2 in the equation of motion (|16[) . Very quickly, in less than an e- folding, nonlinear interactions begin to drive the 
system towards gaussianity. A very simila r behaviou r has been observed in lattice simulations of out-of-equilibrium 
interacting scalar fields during preheating [l65l Il66j | . In the case of rescattering during preheating, the system will 
eventually become gaussian when the fields thermalize. However, in our case the universe is still inflating. As a result, 
nongaussian inflaton fluctuations generated by rescattering are stretched out by the quasi-de Sitter expansion and 
must freeze once their wavelength crosses the Hubble scale. Hence, at late times the PDF docs not become completely 
gaussian, but rather freezes-in with some non-trivial skewness. Within a few e-foldings from the moment of particle 
production the time evolution of the PDF has become completely negligible. 

In order to characterize the nongaussianity of the observable primordial fluctuations, we would like to construct the 
PDF for the curvature perturbation £, including both the contributions from the vacuum fluctuations of the inflaton 
and also from rescattering. To this end, we construct £ using the naive relation ( — —^Scf) and take into account 

both contributions to 8<j> in equation (fT7| . In Fig. [14] we plot the full PDF obtained in this manner, evaluated at very 
late times, well after all relevant modes have crossed the horizon and become frozen. 

Given our numerical results for the PDF of the total observable curvature fluctuation, that is Fig. [TJ] it is straight- 
forward to compute dimcnsionlcss cummulants such as the skewness (£3) and kurtosis (£4) for various values of the 
coupling g 2 . We have summarized our results in Table HV] Note that for g 2 = 0.01 the kurtosis K4 is too small to be 
measured accurately from our HLattice simulations. 

In order to give some sense of the magnitude of the nongaussianity from particle production we have also computed 
an "equivalent /]^£ al " defined by 5k3/(18ct^) where the variance is <7£ = (C 2 ) 1 ^ 2 ~ 10~ 9 / 2 . For a given g 2 this effective 
/]^£ al is the magnitude of fjj l which would be necessary to reproduce the skewness £3 of the IR cascading PDF using 
a local ansatz C = Cg + ISnl [<* - (C 2 )] - 13 

From Table IIVI we see that IR cascading during inflation can generate significant nongaussianity. Even taking 
g 2 = 0.01 (which is compatible with cosmological data for any choice of </>o Q) we still obtain a skewness k 3 — 
—0.006, which is the same value that would be produced by a local model with Jnl^ —53. This "equivalent" local 
nongaussianity is comparable to current observational bounds, and is well within the expected accuracy of future 
missions. This suggests that nongaussian features from particle production during inflation might be observable for 



Our sign conventions for /jvl are consistent with WMAP Q. See [2HI for a discussion of various conventions employed in the literature. 
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FIG. 14: The PDF of the total curvature fluctuation, £, at late times (well after all relevant modes have crossed the horizon 
and frozen). The solid black curve is the exact result from our HLattice simulations and the dotted red curve is a gaussian fit. 
We have also plotted the leading correction to the gaussian result in the Edgeworth expansion, given explicitly by equation 
(24). For illustration, we have chosen g 2 = 0.1 and a standard chaotic inflation potential V(<j>) = m 2 2 /2. 



TABLE IV: Moments of the Probability Distribution Function 



g 2, skewness kurtosis 
(k 3 ) (k 4 ) 


"equivalent" 

rlocal 
JNL 


1 -0.51 0.2 


-4500 


0.1 -0.49 -0.1 


-4300 


0.01 -0.006 < O(10" 3 ) 


-53 



reasonable values of g 2 . 

The "equivalent" /]^£ al values presented in Table IIVI must be interpreted with care. We have included this infor- 
mation only to give a heuristic sense of the magnitude of nongaussianity in our model. It must be stressed that the 
PDF plotted in Fig. Q3] is quite different from the analogous result for local-type nongaussianity. For example, the 
value of the kurtosis (and higher moments) are different, as is the ordering of the cummulants. Moreover, we should 
emphasize that observational bounds on /]^£ al cannot be directly applied to our model since the bispectrum in our 
case is uncorrelated with the vacuum fluctuations and is far from scale invariant. A detailed study of the detectability 
of nongaussianity from particle production will be the subject of a upcoming publication Q . 

Depending on the value of 4>q, the model ((4]) may lead to a variety of observable signatures. As discussed previously, 
0o controls the location of the feature in the primordial power spectrum (|18|) . Nongaussian effects are also localized 
near the same characteristic scale, /cir . If feiR corresponds to scales relevant for CMB experiments, then we predict 
a bump- like feature in the primordial power spectrum, P^(k), and an associated feature in the bispectrum B^(k{) 
with an unusual shape (that will be discussed the next section). A key question is whether the nongaussian feature 
can be observably large in a regime where the power spectrum feature is small enough to be compatible with current 
observations. Preliminary results are encouraging: for g 2 = 0.01 the power spectrum is consistent with all observational 
data [7[ while the skewness of the PDF is rather large. A detailed investigation will require a simple, separable template 
for the bispectrum and will be discussed in a future publication Q . 

On the other hand, we could imagine a scenario where the featu re from IR cascading shows up on smaller scales, 
relevant for Large Scale Structure (LSS) experiments [2ll. Il43l - 145] . In this case our scenario could be probed using 
higher order correlations of LSS probes (such as the galaxy bispectrum) or the abundance of collapsed objects (or 
voids). The latter possibility is interesting since the cluster/void abundance is determined the tails of the PDF and 
may be insensitive to the detailed shape of the bispectrum. Quantitative predictions for observable cluster/void 
abudances require the PDF of the evolved density field, smoothed on some relevant scale 0, rather than the PDF 
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FIG. 15: The shape function S(k,kx2,kxz), defined by ()135p . as a function of the dimensionless quantities £2, £3 which 
parametrize the shape of the triangle. The upper left panel corresponds to k = e _1 fcbum P , the upper right panel is k = fcbum P , 
the lower left panel is k — e +1 fcbum P and the lower right panel is k — e +2 fcbum P - In the IR (k < fcbum P ) the shape of the 
bispectrum is similar to the equilateral shape, however, there is also some support on flattened triangles near k ~ e +1 fcbum P - 
At larger values of k the shape is unlike any other template proposed in the literature. 



of the primordial curvature perturbation (which is plotted in Fig. H~4l . However, we can nevertheless describe the 
qualitative signatures which should be expected. Our model robustly predicts a negative skewness for both the 
curvature perturbation, £, and the density field, Sp/p. Hence, we should expect a decrease in the abundance of the 
largest collapsed objects and an increase in the abundance of the largest voids jl69l Il70j ] . Owing to the localized 
nature of the bispectrum feature, we expect that this effect should show up only when the density field is smoothed 
on a scale close to fciR. 

It is worth mentioning tha t rec ent weak lensing measurement of the dark matter mass of the high-redshift g alaxy 
cluster XMMUJ2235. 3-2557 [l67j | have been construed as a possible hint of nongaussian initial conditions |l6l| . 
Unfortunately, our model does not produce the correct sign of skewness to explain such observations. 



C. Shape of the Bispectrum 

In section IIVI we presented a schematic discussion of how to compute analytically the primordial bispectrum from 
IR cascading, equation ((741 . In Q this expression is evaluated analytically. We find that the bispectrum (1741 from 
particle production has a unique shape that has not been considered in previous literature. In order to describe this 
novel shape we would like to first factor out that strong overall dependence of B c j,(ki) on the size of the triangle. To 
this end we define a "shape function" S(ki) in terms of the bispectrum as follows 

S(fci) = A^- 1 (fc 1 fc 2 fc 3 ) 2 B(fc. i ) (135) 

where B{ki) is related to B^(ki) by equation (|106p and N is a normalization factor to be discussed shortly. The 
function S(ki) has the advantage that the strong k e running of the bispectrum is extracted. Hence, any residual 
scaling behaviour displayed by S(ki ) must be a result of nonlinear interactions. 

Since we expect the 3-point correlation function to be of order P^(k) a natural choice of normalization is N — 

(2tt) 4: V^ where V^ 2 = 5 x 10~ 5 is the usual amplitude of the scale invariant fluctuations from inflation. With 
this choice of normalization our function S(ki) coincides with the quantity Q (kj) / (fci^fcs) which was used to study 
localized nongaussian features from models with steps in the innaton potential in [4ll . |42j . 

Symmetry of the bispectrum under permutations of momenta implies that we can focus only on the region k\ > 
&2 _ fc3, to avoid counting the same configuration twice. Moreover, the triangle inequality implies that 1 — < I 3 -. 
Therefore we can completely specify the shape of the bispectrum for a given size of triangle k by plotting S(k, kx2, kxz) 
in the region X3 < 2; 2 < 1 and 1 — X2 < £3. Because our bispectrum is very far from scale- invariant, it follows that 
this shape function is sensitive to the choice of k. Therefore, in Fig. [15] we choose several representative choices: 
ln(/c/fc bump ) = -1,0,1,2. 

We see that a rich array of shapes are possible: for k < fcbump the bispectrum is qualitatively similar to the 
equilateral model, however, at slightly larger k there is considerable support on flattened triangles also. Note that for 
k > 7.4 fcbump the shape of the bispectrum is extremely unusual and is not easily comparable to any shape that has 
been proposed in previous literature. 
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We find that the bispectrum from IR cascading is, to a first approximation, factorisable in the sense that 
B(ki, fca> ^3) — T[i=iF(ki) where F(k) ~ Cik 3 e~ C2k for some constants C\ and C g ■ T his separable form is im- 
portant since it permits fast algorithms for forecasting, data analysis and simulations |l68j |. 

In this section we have attempted to characterize the nongaussian signature associated with inflationary particle 
production. This signature is rather unique in the literature. Our model predicts uncorrelated localized nongaussian 
features with a unique shape of bispectrum. We have quantified the size of this effect by studying the dimensionless 
cummulants (such as the skewness) and have argued that this nongaussianity can be significant, depending on g 2 . We 
leave a detailed study of the observational constraints on such nongaussianities to future studies [9| . 

D. Multiple Bursts of Particle Production 

Note that our discussion of nongaussianity generalizes easily to the case where there are multiple points 4>i (i = 
0, • • • , n) along the inflaton trajectory where new degrees of freedom \i become massless, such as the model (|107l) . 
Such a construction may be quite natural in the context of brane/axion monodromy inflation [36, 58-60]. In this case 
the nongaussian features in B(ki) from each burst of particle production may superpose to generate a broad-band 
signal. Such a construction may (but need not) be associated with trapped inflation 36]. Depending on the spacing 
of the points 4>i and the couplings g 2 a rich variety of nongaussianities may be possible. In this case that the point 4n 
are sufficiently densely spaced, we expect a nongaussian signal which is close to equilateral but also has some support 
on flattened triangles, consistent with the analysis of [H, [i2|. We leave a detailed discussion to future studies. 

IX. CONCLUSIONS 

In this review we have considered the possibility that some non-inflaton (iso-curvature) particles were produced 
during the observable range of e-foldings of inflation. Inflationary particle production might occur as a result of a 
phase transition, parametric resonance or other non-adiabatic processes. In order to illustrate the basic physics of 
inflationary particle production we restricted our analysis to a simple prototype model with coupling %-(0 — 4>o) 2 X 2 
between the inflaton <p and iso-inflaton x- However, we expect that our qualitative results will apply also to SUSY 
models, gauged interactions, higher spin iso-inflaton fields and phase transitions. 

Models of the type we study have attracted considerable interest recently in connections with trapped inflation, 
trans-Planckian effects and observable features/nongaussianity in the primordial curvature fluctuations from inflation. 
Moreover, such models are quite natural from the microscopic perspective and may be obtained in popular models of 
open string inflation, such as brane/axion monodromy. 

2 

We have shown that inflationary particle production in the model ^-(</> — </>o) 2 X 2 leads to a new mechanism for 
generating cosmological fluctuations. This mechanism is qualitatively different from previous proposals in that we 
do not rely on the quantum vacuum fluctuations of light fields during inflation. Rather, the scenario involves the 
production of massive x particles during inflation, which subsequently rescatter off the slow roll condensate <f>(t) 
to emit bremsstrahhmg radiation of light inflaton fluctuations 6(f). We have studied this dynamics using classical 
lattice field theory simulations, analytical QFT computations and also second order cosmological perturbation theory. 
All of these approaches yield consistent results. We have found that rescattering proceeds with a time scale short 
compared to the expansion time. Moreover, the emission of long- wavelength inflaton fluctuations is very energetically 
inexpensive. The combination of these two effects leads to a rapid build-up of power in IR inflation fluctuations 
shortly after the moment of particle production. This dynamical process is called IR cascading. 

Our numerical and analytical studies of rescattering and IR cascading during inflation may have relevance for 
trapped inflation, preheating, moduli trapping and also non-equilibrium QFT more generally. For instance, we have 
seen that, even with a small number of out-of-equilibrium x particles, multiple rescatterings can nevertheless generate 
long-wavelength 8<f) fluctuations with huge occupations numbers. We have also observed, for the first time, the 
dynamical approach to the scaling regime discovered in [7^, [73| • 

IR cascading during inflation leads to observable features in the primordial cosmological fluctuations. In particular, 
this process generates a bump-like contribution to the primordial scalar power spectrum. This signature is very 
different from what would be obtained in a model with transient violation of slow roll during inflation (such as a 
step-like feature in V{<p)), contrary to some claims in the literature. 

We have studied the observational constraints on bump-like features from inflationary particle production during 
inflation. We found that relatively large distortions, of order 10% of the usual scale invariant vacuum fluctuations, 
are compatible with current data. We have derived observational bounds on the coupling g 2 for a given (fto, which 
play a crucial role in determining the detectability of nongaussianity from particle production. Our observational 
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bounds on particle production during inflation have implications for brane/axion monodromy inflation models and 
other microscopic constructions. 

IR cascading also has a nontrivial impact on nongaussian statistics, such as the bispectrum. The model 4j- (4>— 4>o) 2 X 2 
leads to a very novel nongaussian signature: uncorrelated, localized nongaussian features with a unique shape of 
bispectrum. For reasonable values of the coupling, g 2 < 0.01, this new kind of nongaussianity may be detectable in 
future missions. 

The nongaussian signature predicted by inflationary particle production is rather unusual, as compared to other 
models of inflation which are frequently studied in the literature. However, the underlying field theory description is 
extremely simple and rather generic from the low-energy perspective. In order to obtain observable nongaussianity it 
was not necessary to fine-tune the inflaton trajectory or appeal to re-summation of an infinite series of high dimension 
operators. Indeed, the only "tuning" which is required for our signal to be observable is the requirement that <fi — 0o 
during the observable range of e-foldings. 

There are a variety of directions for future studies. From the theoretical perspective, it would be interesting 
to explicitly generalize our results to more complicated models with particle production during inflation (such as 
SUSY models, fermion \ fields and phase transitions). There are also a wide range of interesting phenomcnological 
possibilities. Varying the location of the feature we can have a variety of possible signatures for the CMB and LSS. 
We expect that IR cascading will also have implications for the spectrum of gravity waves from inflation and also 
primordial black holes. We can also imagine superposing multiple bursts of particle production to obtain an even 
richer variety of signatures. We leave these possibilities for future investigation. 
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APPENDIX A: Backreaction Effects 



Once <f> rolls past the massless point (f> — (f>o a, gas of x particles is produced with occupation number given by (|13[) . 
This gas costs energy, which must be drained from the condensate 4>(i). Hence, in order for the system to conserved 
energy, the inflaton must slow down slightly. As discussed in [(| Ea . l46l loH |68| , this slowing down can be studied 
using the mean field equation (fT4")l . Using the solutions (|47l) to compute the vacuum average (x 2 ) we find 

(2ttY a J (t) \4> - 0o I 

The step function Q(t) reflects the fact that the impact of particle production is felt only after <j> passes through 0o 
and a -3 comes from the volume dilution of non-relativistic particles. The final term in (|A-1| may be interpreted as 
a quantum correction to the effective force due to particle production. (Note that we have implicitly subtracted the 
Coleman- Weinberg potential, which may be justified by the assumption of softly broken SUSY.) 

The solutions of (|A-1I) display the expected slowing-down behaviour and have been studied analytically in [H, |H, 
l6l| . In Q this slowing-down was studied using inhomogeneous lattice field theory simulations and the results were 
found to be compatible with the mean field treatment (see Fig. [5]). Here, we consider simple energetic arguments 
in order to clear up some common misconceptions. The effective inflaton potential, including the effects of particle 
production, is 



kl \<t>-M 
' (2tt) 3 



V eS = V{<t>)+g@{t) * |y 3 yu| (A-2) 



For t > K 1 the X particles are non-relativistic and their energy density is dominated by potential (rather than kinetic) 
energy. Hence, we have 
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with TV = Ht the number of e-foldings measured from the moment when (ft = <j)Q. Shortly after particle production 
this energy density grows, corresponding to the fact that as <f> moves away from (f> the x particles become ever more 
massive. However, this growth in the energy density cannot continue forever. At N = 1/3 the energy density p x peaks 
and at later times it decays exponentially as e~ 3N , corresponding to the volume dilution of the massive x- Thus, the 
energy density in the produced x is always bounded as 



Px < 



1 



247r 3 e H 



(A-4) 



On the short time scales relevant for particle production, | AN\ < 1, it is still sensible to talk about energy conservation. 
The energy in produced x particles must therefore be balanced by a dip in the kinetic energy of the inflaton. (Intuitively 
this is to be expected since the "extra" term in (|A-1I) represents a force which tends to pull the inflaton back towards 
the point <fi = 4>q.) To get a sense of the magnitude of this effect, let us compare p x to the initial kinetic energy of 
the inflaton 



Using (fTTj) it is easy to see that 



Kb 



2ff 2 



^- < 0.35 g b ' 2 « 1 
K- 



(A-5) 



(A-6) 



Hence, the total energy density that goes into x particles is small compared to the inflaton kinetic energy, and hence 
the velocity dip must also be a small effect, roughly &-4>/4> < 0.18 g 5 ^ 2 - This simple estimate is consistent with a more 
quantitative treatment [f|. 

Note that, since inflation is driven by potential energy, we have V 3> K m 3> p x . Hence, even during particle pro- 



duction the expansion rate H will still be dominated by the classical inflaton potential: H = JV/ (3Mp). Therefore, 

a single burst of particle production will not terminate inflation. Nor will this burst de-cohere the condensate since 
the occupation number of produced particles (fl3| is always less than unity. (We again remind the reader that, unlike 
the case of broad band resonant preheating after inflation, we have only a single burst of particle production.) 

Note that throughout this appendix we are assuming standard slow roll inflation. In particular, our argument does 
not apply for trapped inflation [36| where the effect of backreaction on H and <fi(t) is much more significant. 



APPENDIX B: Detailed Computation of P(k) 

In this appendix we discuss some of the technical details associated with the computation of the renormalized power 
spectrum (jTTj) . First, notice that using (|43p and (|4"T)) we can write the quantity appearing in each renormalized Wick 
contraction as 



1 1 1 



Tlk COS 



klt\r) klt 2 {r>) 



+ \/nk\/^ + n k sir 



k 2 t 2 (r) *#V) 



where the occupation number is defined by (fl"3"]) . Plugging (|B-1[) into (JTTJ) we find 



(B-l) 



g 2 k 3 



d krik-k'n-k' x / dr dr 
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(B-2) 



Notice that the time and phase space integrations in (|B-2|) decouple. This is the key simplification which makes an 
analytical evaluation of this expression tractable. Let us consider these integrations separately. 



4G 



B.l. Time Integrals 



All of the integrals over conformal time that appear in (|B-2|) can be expressed in terms of two characteristic integrals 
which we call 1\ and 7 2 - Explicitly, these are defined as 



h(k,r) = t~t" f dr'G k {T-r')e lk " t2 ^ 
a (V J 

h(k,T) = / dT'G k (T-T') 



(B-3) 
(B-4) 



The second characteristic integral, / 2 , can be evaluated analytically. However, the resulting expression is not partic- 
ularly enlightening. Evaluation of I±, on the other hand, requires numerical methods. 

Let us now show how the various integrals appearing in (|B-2| may be re- written in terms of 1%, / 2 . First, consider 
the first line of (|B-2j) where the following integral appears: 



a{r) a{T) 
\h(k,r)\ 2 , h{Krf 



klt\r<) klt\r") 



2 2 

Next, consider the second line of (|B-2j) where the following integral appears: 

tklt\r>) , fc 2 t 2 (r") 
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Finally, consider the fourth line of (|B-2j) where the following integral appears 

"A^V) fc 2 i 2 (r") 
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In the expressions (|B-6p and (|B-7[) the notations Re and Im denote the real and imaginary parts, respectively. 



B.2. Phase Space Integrals 



As a warm-up to the subsequent calculation consider the following integral: 



d 3 k'n a k „ k ,n b k , 



d 3 k'exp [-a7r|k-k'| 2 //c 2 ] exp [-67r|k'| 2 /A; 2 ] 



(a + hf/ 2 



exp 



ab nk 2 
i + b ki 



(B- 



This formula is valid when a, b are positive real numbers. Notice that this expression is symmetric under interchange 
of a and b. 

The phase space integral in the first line of (|B-2j) is computed by a trivial application of the identity (|B-8|) : 



d 3 k'n k ^ k ,n k i 



= -7rfc 2 /(2fc?) 



2V2 



(B-9) 



However, the remaining phase space integrals appearing in (|B-2[) cannot be obtained exactly in closed form because 
they contain terms like \/l + n k > where the gaussian factors appear under the square root. In order to deal with such 
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expressions, we note because n k <C 1 over most of the domain of integration, it is reasonable to replace \/l + n&' 
1 + nt~i /2. Let us now proceed in this manner. The phase space integral on the second line of (|B-2I) is: 



^k'y/nk-k'rik'^l + rik-k'Vl + n k > 
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Finally, consider the phase space integral on the third line of (IB-2|) 
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We have verified the formulae (|B-10IB-11[) numerically. In both cases that the numerical results agree with these 
semi-analytical expressions up to the percent level. 

We can now, finally, insert the results (|B-5IB-6IB-7[) and (|B-9IB-10IB-11[) into the expression (|B-2[) . Doing so, we 
arrive at our main analytical result, which is equation (|71l) . 
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